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Abstract 


This  report  summarizes  aeroservoelastic  tailoring  studies  in  which  adaptive  material 
actuators  are  used  to  control  structural  deflection  of  aeroelastic  systems.  The  actuator 
must  furnish  enough  control  of  a  system  to  make  the  control  of  the  phenomenon  feasible. 
Specific  research  problems  considered  are:  choice  of  the  actuator  material  for  effective 
control;  geometric  arrangement  of  actuators  for  active  control;  and  optimum  coverage  of 
surface  panels  for  effective  control;  a  controller  design  process  to  improve  efficiency  of 
the  selection  process.  A  finite  element  method  is  developed  to  model  actuator  and  sensor 
output  for  plate-like  actuators  and  its  use  is  illustrated  for  wing-like  configurations  to 
demonstrate  the  benefits  of  orthotropic  material  actuators.  The  problem  of  optimum 
actuators  to  supply  deflection  of  panels  for  wing  surfaces  is  examined  to  determine 
optimality  criteria  for  such  panels  and  to  use  strain  energy  as  a  guide  for  efficient  use  of 
actuator/host  plate  combinations.  Finally,  a  method  of  controller  design  is  suggested  to 
determine  the  limits  of  actuator  control  effectiveness  for  a  typical  section  whose  response 
to  random  atmospheric  turbulence  is  to  be  minimized. 
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Section  1.0  •  Introduction 


1.1  Purpose  and  Scope 

This  report  summarizes  studies  conducted  between  1  October  1991  and  30  September 
1993.  The  purpose  of  this  research  was  to  develop  fundamental  understanding  of  the 
features  of  successful  integrated  aeroservoelastic  design  when  active  materials  are  used  as 
actuators  to  control  the  deformation  of  lifting  surfaces.  Key  features  of  this  study  focus 
on  the  development  of  appropriate  models  and  identification  of  effective  analytical 
techniques. 

The  first  problem  studied  is  the  process  of  actuator  material  selection  and 
development  new  a  finite  element  model  and  optimization  studies  to  determine  the 
conditions  for  optimality  of  actuator  host  plate  configurations.  This  latter  effort  led  to  a 
strain  energy  approach  for  actuator  design.  Secondly,  the  research  examined  the 
conditions  for  the  effective  integration  of  control  estimation  theory  into  the  aeroelastic 
process  when  active  materials  were  used.  This  effort  produced  results  summarized  in 
Section  3.0. 

1.2  Summary  of  Research  Efforts  -  Publications  and  Education 

As  the  result  of  this  research  funding,  one  Master’s  Degree  student,  Ms.  Tamara 
Leeks,  was  funded  for  her  early  research  work,  as  a  student  of  the  Principal  Investigator, 
Professor  Weisshaar.  She  later  received  a  National  Science  Foundation  Fellowship  to 
continue  her  studies  into  the  Ph.D.  program.  Ms.  Leeks  received  her  M.S.  degree  in 
December  1993  and  is  now  a  Ph.D.  student  at  Purdue,  continuing  her  research  on 
piezoelectric  actuators.  Ms.  Leeks  and  Professor  Weisshaar  will  present  a  paper  at  an 
AIAA  Conference  in  April  1993  based  on  results  obtained  under  this  AFOSR  grant 
sponsorship. 

Professor  Mario  Rotea,  Faculty  Associate,  also  furthered  his  research  efforts  in  the 
design  of  effective  controllers  for  active  aeroelastic  control.  He  is  the  author  of  the  work 
covered  in  Section  3.  He  is  now  a  National  Science  Foundation  Presidential  Young 
Investigator. 

This  funding  also  supported  Professor  Weisshaar’s  efforts  to  develop  advanced  finite 
element  analysis  models  by  allowing  leveraging  funding  for  collaboration  with  a 
Fulbright  visiting  scholar.  Professor  M.N.  Abdul-Wahed.  The  results  of  this 
collaboration  are  contained  in  the  Appendix  as  a  paper  presented  at  a  national  conference. 
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1.3  Controlled  structure  synthesis 

Lifting  surfaces  such  as  wings  and  stabilizer  surfaces  are  flexible  and  subject  to 
dynamic  response  phenomena  including  flutter,  divergence,  gust  response  and  buffeting. 
Because  of  this,  hydraulically  driven  aerodynamic  control  surfaces  such  as  ailerons,  are 
used  to  generate  aerodynamic  loads.  These  controlled  loads  respond  to  the  sensed  motion 
of  the  parent  surface  to  cancel  out  objectionable  response  such  as  high  stresses  or 
structural  instabilities. 

Lifting  surface  response  in  flight  is  complicated  by  complex  interactions  between  the 
flow  field  around  the  wing  and  the  response  of  the  wing  itself.  This  interaction  is  called 
aenoelasticity.  The  use  of  servo  controllers  to  control  aeroelasdc  response  is  called 
aeroservoelasticity.  Two  types  of  aero servoe las ticity  problems  were  addressed  in  this 
study.  The  first  problem  is  the  efficiency  of  piezoelectric  actuators  to  create  and  shape 
the  lifting  surface  loads  and  structural  response.  These  actuators  are  imbedded  into  the 
structural  surface  so  that  they  generate  load  and  carry  stresses.  The  second  problem 
development  is  the  modeling  of  the  host  structure  and  its  actuators  so  that  effective 
designs  can  be  evaluated. 

The  goal  of  adding  an  active  control  system  to  a  structure  is  to  improve  performance 
of  the  original  system  by  the  addition  of  sensors,  control  processors  and  actuators.  An 
even  more  important  goal  is  to  achieve  performance  that  cannot  be  done  in  any  other 
practical  way.  For  instance,  there  may  be  no  uncontrolled  structure  with  the  requirements 
placed  on  it  that  fulfills  the  design  goals.  A  central  issue  in  the  design  of  a  controlled 
structure  is  the  selection  and  placement  of  actuators  and  sensors  on  and  within  the  system 
to  be  controlled. 

A  well-designed  controlled  structure  is  tolerant  of  imperfections  in  the  system.  This 
feature  is  called  robustness.  The  controlled  system  should  also  have  good  performance  in 
rejecting  system  disturbances  (gusts  for  instance).  This  means  that  it  should  attenuate  the 
effects  of  these  disturbances.  In  the  controls  community  this  attenuation  is  called 
regulation. 

When  designing  a  control  system,  care  must  be  taken  to  use  a  technique  that  monitors 
the  signals  sent  by  the  controller  computer  to  the  actuator  so  that  they  are  not  too  large  so 
that  the  actuator  limits  are  exceeded.  In  the  case  of  an  aileron  controller,  this  means  that 
the  angle  of  aileron  rotation  should  not  be  too  large.  In  the  case  of  a  piezoelectric  device, 
it  means  that  the  electric  field  applied  to  the  active  material  should  not  be  too  large. 

One  major  problem  with  control  law  evaluation  and  relating  control  law  requirements 
to  actuator  choices  is  that  structural  control  research  has  been  restricted  to  a  control  law 
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that  receives  signals  from  the  sensors  and  changes  them  into  signals  to  the  actuators  that 
are  in  a  specific  form.  Restricting  the  form  of  a  control  law  or  confining  the  study  to  only 
one  type  of  solution  algorithm  is  undesirable  for  two  reasons.  First  of  all,  when  the 
designer  is  done,  he  has  not  answered  the  central  question  that  the  structural  engineer  is 
interested  in.  This  question  is  Have  I  designed  the  best  possible  combination  of 
structure  and  actuator  (including  size  and  location  on  the  structure)  so  that  the  result  does 
the  job  efficiently?"  Secondly,  when  the  process  is  complete  it  is  not  certain  that  the 
controlled  structure  represents  the  best  combination  of  host  structure  and  sensor  actuator 
design. 

The  method  examined  in  this  study  was  not  restricted  to  a  single  method  of  control 
law  design.  In  fact  it  will  reveal  only  what  the  best  configuration  is,  but  not  the  necessary 
control  law.  Said  differently,  this  new  method  will  tell  us  what  the  limits  of  performance 
are,  but  not  how  to  formulate  the  control  law  to  get  there.  This  is  a  key  feature  of  this 
method.  This  type  of  information  is  more  valuable  to  the  structural  engineer  who  is  not 
concerned  with  the  control  law  algorithm  or  its  development 

1.4  Issues  with  integrated  active  material  sensor/actuators  for  aeroservodastidty 

Adaptive  structures  use  actuators  to  create  changes  in  structural  design  shape  or 
damping  to  respond  to  changing  mission  needs  and  performance  requirements.  Actuators 
may  be  articulated  mechanical  devices  such  as  ailerons  and  flaps  or  they  may  be  "smart" 
materials  embedded  in  the  structure.  In  the  case  of  atmospheric  flight  vehicles,  active 
material  actuators  embedded  in  the  structure  can  control  structural  shape  and 
aerodynamic  loads. 

Using  sensors  and  feedback  control  algorithms,  these  actuators  also  can  change  wing 
stiffness  and  control  the  stability  of  an  entire  lifting  surface  to  maneuver  and  trim  an 
aircraft,  reduce  gust  loads  or  enhance  the  stability  of  the  structure.  In  some  cases,  this 
can  be  done  with  less  weight  and  at  reduced  overall  cost  compared  to  conventional 
structures.  At  present,  these  concepts  rely  on  aeroelasticity,  that  is,  the  mutual  interaction 
between  aerodynamic  loads  and  structural  deflections,  to  create  favorable  active  load 
control. 

The  design  of  controlled  structures  involves  trade-offs  such  as  structural  stiffnesses, 
actuator  choices  and  the  location  of  actuators.  In  aeroservoelastic  design,  this  selection 
process  determines  the  trade-offs  between  conventional  control  surfaces  and  active 
materials.  Which  is  more  effective  in  achieving  the  design  goals,  the  active  material  or 
the  aerodynamic  surface? 
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Aerodynamic  loads,  and  the  local  pressures  that  create  these  loads,  depend  on  the 
surface  shape  on  which  they  act  Surface  panels  may  be  flat  or  curved  and  are  designed 
to  provide  aerodynamic  shape  and  to  guarantee  structural  integrity.  Re-shaping  smooth 
aero/structural  surfaces  to  change  the  pressure  distribution  is  done  by  bonding  or 
otherwise  attaching  thin  actuators  to  the  inside  surfaces  of  structural  panels  to  create  an 
asymmetric  configuration  that  will  bend  on  command. 

Thin  plate-like  or  lattice  reinforced  panels  with  embedded  self- straining  actuators 
such  as  shape  memory  materials  or  piezoelectric  materials  have  been  proposed  for 
aerodynamic  control  concepts  that  include  actively  controlled  panels  to  reduce  transonic 
drag  and  active  panel  elements  to  increase  supersonic  panel  flutter  speed.  For  transonic 
drag  reduction,  the  deformation  of  a  panel  on  the  upper  surface  of  a  supercritical  airfoil 
can  change  the  flow  field  and  shock  wave  intensity  to  reduce  drag  on  command. 

Panel  flutter  suppression  v;th  piezoelectric  actuators  and  shape  memory  alloy 
actuators  is  unique  in  that  no  articulated  device  exists  to  do  the  same  task.  In  supersonic 
flow,  dynamic  oscillations  can  be  reduced  by  placing  thin  actuators  on  the  panel  surfaces 
to  change  the  frequencies  of  the  panels  on  command. 

One  problem  with  active  panel  concepts  is  the  difficulty  finding  a  design  combination 
to  give  large  enough  panel  out-of-plane  deformations  to  create  the  required  changes  in 
aerodynamic  forces.  Without  deflections  of  the  order  of  a  panel  thickness  (or  even  more), 
controlling  the  size  and  position  of  the  aerodynamic  forces  is  marginal. 

A  desirable  actuator,  such  as  one  using  today's  piezoelectric  materials,  can  not  create 
significant  bending  deformation  of  panels  unless  the  host  panel/actuator  combination  is 
tailored  to  extract  every  bit  of  electro- mechanical  efficiency  out  of  the  configuration.  An 
emphasis  on  efficiency  naturally  leads  to  considerations  of  formal  optimization  that 
includes  a  design  objective  and  design  variables.  However,  before  formal  optimization 
can  proceed,  we  must  select  our  design  variables  and  determine  the  sensitivity  of  the 
design  objective  to  these  design  variables. 

Section  2  of  this  report  reviews  pre-optimization  studies  to  explore  the  interaction 
between  actuator  self-straining  ability  and  bending  stiffness,  thickness  and  planform 
coverage  and  the  host  panel  bending  stiffness  and  aspect  ratio.  The  purpose  of  the 
actuator  is  to  produce  large  bending  deflection.  The  intent  of  the  study  is  to  identify 
effective  panel/actuator  combinations  and  understand  why  some  combinations  are  more 
effective  than  others. 
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Section  2.0  -  Actuator  optimization  and  analysis  studies 
2.1  Purpose  and  scope  •  aeroeiastic  load  control 

This  section  reviews  the  purpose  and  results  of  to  studies  related  to  actuator/host 
structural  matching  and  optimization  for  performance  of  the  actuator.  Control  of 
aerodynamic  loads  and  structural  response  of  wings  and  other  lifting  surfaces  is  not  a  new 
concept  Proposals  to  control  wing  response,  in  particular  gust  response  are  over  thirty 
years  old.  Ride  quality  enhancement  was  used  on  the  B-70  Valkerie  Supersonic  Bomber. 
The  same  concept  is  used  today  on  the  B- 1  bomber,  to  allow  it  to  conduct  low  level ,  high 
speed  penetration  missions.  The  Lockheed  L-1011  uses  active  inboard  flaps  to  reduce  the 
seventy  of  turbulence  on  wing  fatigue  life  and  on  the  ride  quality  felt  by  passengers. 

Suppression  of  wing  flutter,  a  dynamic  structural  instability  created  by  unfavorable 
interaction  between  unsteady  aerodynamic  loads  and  structural  vibrations,  was 
demonstrated  in  the  mid-1970's.  This  flutter  suppression  was  possible  by  using  feedback 
control  to  drive  aerodynamic  surfaces  such  as  ailerons  and  leading  edge  slats  to  damp  out 
motion. 

A  piezoelectric  material  is  a  material  which,  when  subjected  to  mechanical  load, 
develops  an  electric  charge  proportional  to  the  resulting  mechanical  stress.  Conversely, 
this  material  will  deform  or  strain  when  an  electric  field  is  applied.  This  deformation,  in 
particular  positive  or  negative  strain,  depends  on  the  polarity  of  the  applied  field.  The 
electric  field  is  generated  by  imposing  a  voltage  across  the  material  so  that  the  field  is 
proportional  to  the  applied  voltage  divided  by  the  distance  between  the  electrodes. 

Active  materials,  in  particular  piezoelectric  material  actuators,  have  been  proposed  as 
controllers  to  accomplish  much  the  same  mission  as  aerodynamic  surfaces.  How 
successful  and  in  what  areas  these  material  actuators  will  be  able  to  replace  or  augment 
aerodynamic  controllers  such  as  ailerons  is  still  very  much  in  doubt.  There  are 
advantages  to  creating  so-called  "solid  state  wings"  that  have  no  external  hinges  or 
mechanical  devices.  For  military  applications  there  is  the  advantage  o*  stealth.  Any  gaps 
or  openings  that  appear  in  the  wing  tends  to  increase  size  of  the  radar  return  from  an 
aircraft  Another  advantage  is  that  the  control  is  distributed  over  the  entire  wing,  making 
control  of  some  aeroeiastic  phenomena  more  precise.  Lastly,  there  are  some  aeroeiastic 
phenomena,  such  as  panel  flutter,  for  which  no  mechanical  device  is  feasible. 
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2.2  Summary  of  finite  element  development  and  studies 


The  use  of  actuators  to  control  aeroelastic  deformation  also  depends  on  their  ability  to 
decouple  two  characteristic  types  of  deformation,  bending  and  twist.  This  modal 
decoupling  is  very  important  to  aeroelastic  load  control  because  both  types  of 
deformation  create  aerodynamic  loads,  but  we  may  not  always  want  these  loads  to  occur 
simultaneously. 

The  study  of  aeroelastic  decoupling  requires  an  accurate  analytical  model,  one  that 
can  model  arrays  of  anisotropic  actuators  oriented  at  angles  skewed  with  respect  to 
conventional  structural  axes.  It  also  requires  realistic  modeling  of  an  aircraft  wing 
structure.  This  model  may  be  a  plate-like  configuration  or  the  closed  box-beam 
arrangement  common  to  all  efficient  aeronautical  structures. 

Because  of  the  interest  in  effective  actuator  combinations  and  realistic  structures,  a 
finite  element  method  was  developed  to  model  potential  structure/actuator  configurations. 
The  finite  element  model  consists  of  three-dimensional  isoparametric  solid  elements 
(bricks)  that  allow  modeling  of  tailored  piezoelectrics  with  skewed  actuator/sensor  axes. 
This  scheme  also  allows  the  representation  of  an  anisotropic  host  structure  and  can 
account  for  material  and  stacking  geometry  through  the  element  thickness.  Using  this 
finite  element  representation,  it  is  shown  that  anisotropic  piezoelectric  actuators  can 
create  sufficient  twisting  and  bending  to  control  aerodynamic  loads  on  a  wing. 

To  illustrate  the  application  of  the  finite  element  model  developed  in  the  course  of 
this  research,  this  method  was  used  to  examine  four  different  actuator/host  plate 
combinations  to  illustrate  orthotropic  actuator  control  and  structural  response.  In  one 
case  the  deflection  of  a  bimorph  plate  constructed  of  two  PVDF  layers  so  that  the  top 
layer  acts  as  a  distributed  actuator  while  the  bottom  layer  acts  as  a  sensor  was  examined. 
In  another  case,  a  rectangular  aluminum  plate  was  sandwiched  between  two  piezoceramic 
layers  and  its  deflection  analyzed  to  see  if  bending  and  torsional  deflection  could  be 
separated.  A  similar  case  considered  the  same  aluminum  plate,  but  with  two  layers  of 
tailored  piezopolymers  to  assess  the  effects  of  piezoelectric  anisotropy.  Finally,  a  steel 
plate  was  sandwiched  between  two  PV  JF  elements  with  off-axis,  mirror  symmetry  to  test 
the  ability  of  a  thick  layer  piezoelectric  actuator  to  deform  the  plate. 

The  result  of  these  studies  was  a  paper  presented  at  the  North  American  Conference 
on  Smart  Materials  and  Structures,  held  in  Albuquerque,  N.M.  in  February  1993.  This 
paper  is  included  in  the  appendix  to  this  report  This  paper  concluded  that:  , 

•  Electromechanical  anisotropy  is  an  effective  and  highly  desirable,  if  not  essential,  to 
aeroelastic  control.  Efforts  should  be  directed  towards  developing  more  effective 


6 


anisotropic  materials  with  larger  (and  unequal)  strain  or  stress  constants  (d3i  or  e32) 
and  a  larger  modulus  of  elasticity. 

•  The  finite  element  method  can  provide  an  accurate  solution  for  thin  and  thin-walled 
structures  microactuated  by  thin  piezoelectric  elements.  This  method  can  readily 
account  for  the  presence  of  webs  and  stiffeners,  either  active  or  passive,  to  study  new 
configurations. 

Future  work  will  concentrate  on  using  this  finite  element  model  to  assess  aeroelastic 
effects  on  thin  plates  in  supersonic  flow  and  built-up  wings  in  transonic  flow.  Such 
configurations  can  be  used  to  suppress  flutter,  control  shock  wave  formation  and  to 
reduce  drag.  In  addition,  this  finite  element  formulation  can  be  used  for  optimization 
studies  to  determine  effective  actuator  geometry  and  locations  and  to  help  to  understand 
issues  related  to  integrating  these  actuators  with  host  structures. 

2.3  Panel  actuator  combinations  for  maximum  deflection  -  optimality  criteria 

Another  study  conducted  with  this  funding  examined  the  features  of  a  self¬ 
straining  actuator  mounted  on  one  side  of  a  flat  panel  to  cause  the  largest  deflection.  A 
Rayleigh-Ritz  model  was  developed  to  compute  inplane  and  bending  deflections  of  a 
plate  with  an  actuator  covering  only  part  of  the  area. 

It  was  shown  that  it  is  important  to  include  the  coupling  terms  and  the  in-plane 
actuator  forces  to  accurately  model  the  problem.  A  Rayleigh-Ritz  trigonometric  series 
model  was  developed  to  analyze  the  problem  and  to  solve  for  deflections  and  stresses. 
This  model  allows  the  actuator  to  be  placed  at  any  location  on  the  panel.  It  was  found 
that  the  method  gave  us  excellent  results,  very  close  to  NASTRAN  finite  element 
analysis. 

For  aluminum  panels  with  aspect  ratios  between  1.0  and  1.5,  the  best  actuator  has  a 
thickness  ratio  of  0.6  and  covers  65%  of  the  panel  area.  For  a  panel  with  an  aspect  ratio 
of  3.0,  the  best  actuator  has  a  thickness  ratio  of  0.6  and  covers  70%  of  the  area.  For  the 
panels  with  aspect  ratios  of  1.0  and  1.5,  the  actuator  with  the  same  aspect  ratio  as  the 
panel  produced  the  largest  deflection.  For  the  panel  of  aspect  ratio  3.0,  an  actuator  with 
an  aspect  ratio  slightly  smaller  than  that  of  the  panel  produced  the  most  deflection. 

For  a  plate  with  cylindrical  bending,  when  the  strain  energy  is  plotted  against 
thickness  ratio,  an  inflection  point  occurs  at  the  point  where  the  actuator  thickness  creates 
the  most  center  panel  bending  deflection.  When  the  panel  is  simply  supported  all  around, 
the  plot  of  strain  energy  vs.  actuator  thickness  ratio  becomes  nearly  linear  at  the  point 
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corresponding  to  the  most  effective  actuator  thickness.  This  indicated  that  increasing  the 
thickness  of  the  actuator  was  beneficial  until  the  slope  of  the  strain  energy  curve  reached 
its  smallest  value.  Adding  more  thickness  to  the  actuator  beyond  this  point  increases 
stiffness  more  than  it  increased  the  applied  moments,  and  does  not  increase  the  deflection 
that  is  obtained  in  the  panel. 

More  precise  results  for  the  actuator  characteristics  could  be  calculated  if  an 
optimization  scheme  were  used  in  conjunction  with  the  model  already  developed.  Since 
only  a  finite  number  of  combinations  were  tried  in  this  study,  the  best  actuator  was  found 
within  the  limits  of  the  study.  Also,  the  optimal  actuator  could  be  found  by  taking  into 
consideration  the  weight  that  is  added  for  a  larger  actuator.  A  larger,  thicker  actuator 
does  not  always  produce  significantly  more  deflection  than  a  smaller,  lighter  actuator. 

This  study  was  also  limited  to  rectangular  actuators  placed  at  any  location  on  the 
panel.  Actuators  of  shapes  other  than  rectangular  should  also  be  considered.  Rectangular 
actuators  produce  high  stresses  at  the  edges  and  comers.  The  comer  stresses  might  be 
avoided  if  the  actuator  shape  were  changed.  The  deflection  produced  by  elliptical 
actuators  or  other  actuator  shapes  should  be  examined, 

2.4  Summary 

The  analysis  necessary  to  confidently  produce  optimized  actuators  for  aerodynamic 
and  aeroelastic  control  has  been  improved  because  to  the  results  of  the  studies  just 
reviewed.  Technical  papers  that  detail  this  research  are  contained  in  the  Appendix  to  this 
report.  This  research  area  will  also  be  expanded  to  be  the  topic  of  Ms.  Leeks  Ph.D. 
research. 


3.0  -  Active  control  of  a  typical  section 

This  following  section  summarizes  the  analysis  and  methodology  developed  to  use  to 
evaluate  the  effectiveness  of  aeroelastic  adaptive  material  actuators.  The  example  chosen 
is  the  typical  section.  The  so-called  typical  section,  as  originally  developed,  represented 
an  airfoil  section  on  an  unswept  wing  whose  dynamic  behavin’  was  typical  of  the 
response  of  the  entire  wing.  This  wing  section  is  shown  in  Figure  2  on  the  next  page. 
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The  essential  structural  features  of  this  wing  model  are  its  two  uncoupled  natural  frequen¬ 
cies,  one  a  torsion  reference  frequency,  u>a  -  the  other  a  bending  or  plunge  reference  frequency, 
-V  These  frequencies  are  due  to  elastic  restraint  provided  by  the  idealized  springs  depicted 
in  Figure  2.  An  airstream  impinges  on  the  airfoil  and  generates  lift  by  vortex  action.  The 
amount  of  lift  is  proportional  to  the  incidence  or  angle  of  attack  of  the  surface.  This  lift  is 
regarded  as  being  generated  instantaneously,  with  no  delays  or  lags  between  the  aerodynamic 
forces  and  the  displacements. 


Figure  2:  Two-dimensional  aeroservoelastic  model  [2], 

The  airfoil  mass,  its  moment  of  inertia  and  its  position  relative  to  the  so-called  shear 
center  of  the  section  determine  the  values  of  the  actual  airfoil  natural  frequencies,  which 
are  close  to  u>„  and  w*  when  the  airstream  is  not  present.  These  frequencies  result  from 
dynamic  coupling  between  the  two  uncoupled  section  frequencies.  The  aerodynamic  force 
and  moment  are  referenced  to  a  point  of  action  called  the  aerodynamic  center,  located  at 
the  1/4  chord  aft  of  the  leading  edge,  as  shown  in  Figure  2. 

Active  control  of  this  section  is  furnished  by  some  combination  of  three  independent 
actuators.  The  first  of  these  controls  is  the  conventional  aileron  indicated  in  Figure  2.  This 
aileron  has  limited  effectiveness  because  it  will  cause  the  airfoil  to  stall  and  lose  lift  if  it 
rotates  too  far  with  respect  to  the  airfoil.  This  rotation  is  represented  as  0  in  Figure  2  and 
its  limit  is  0O. 

The  second  type  of  control  comes  from  an  adaptive  piezoelectric  material  that  applies 
an  internal  force  to  the  airfoil  at  the  shear  center  and  also  creates  an  active  torsional  mo¬ 
ment  proportional  to  torsional  displacement.  This  control  is  furnished  by  adaptive  material 
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actuators  imbedded  m  the  cover  skins  and  the  spars  of  the  wing. 

The  equations  of  motion  for  the  airfoil  section  are  written  in  nondimensional  form  as 
follows: 
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Equation  1  includes  the  aileron  controller,  the  inertia  matrix  M  and  the  stiffness  matrix 
Ki,  as  discussed  previously.  It  also  includes  the  matrix  Kj  that  accounts  for  aerodynamic 
forces  that  are  functions  of  airspeed  U  the  displacements,  h{t)  and  a(t).  The  plunge  dis¬ 
placement  is  divided  by  the  semi-chord  dimension  b  =■  c/2  to  obtain  the  nondimensional 
signal  A(t)  =  h(t)/b.  Equation  1  contains  aerodynamic  section  coefficients  appropriate  to 
the  wing. 

Flow  unsteadiness,  an  essential  feature  of  models  for  future  study,  can  be  added  by 
modifying  the  Kj  matrix  to  include  additional  degrees  of  freedom  or  lag  states.  However, 
inclusion  of  such  terms  now  will  unnecessarily  complicate  the  problem. 

The  presence  of  adaptive  materials  is  found  in  the  terms  k  and  g  in  Eqn.  1.  These  terms 
are  found  on  the  right  hand  side  of  that  equation.  The  signal  k(t )  models  the  ability  of  active 
materials  to  control  the  bending  stiffness  of  the  section.  This  signal  is  bounded  and  may 
not  exceed  the  value  k 0  because  of  material  limitations.  Similarly,  the  signal  g(t)  models  the 
ability  of  the  active  materials  to  control  the  torsional  stiffness  of  the  wing.  This  signal  can 
be  either  independent  of  k(t)  or  related  to  k(t).  Here,  we  assume  that  g(t)  does  not  depend 
on  k(t)  and  that  g{t)  is  also  bounded  and  cannot  exceed  the  value  ge.  Both  k0  and  ge  are 
functions  of  actuator  placement,  material  properties  and  geometry. 

The  equation  of  motion  also  contains  the  influence  of  a  random  gust  field.  This  influence 
is  represented  by  the  term  wt,  representing  the  usual  Dryden  gust  field  used  in  aeroelastic 
design.  The  signal  w,(t )  is  the  vertical  velocity  of  the  gust  and  is  generated  by  the  following 
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second  order  dynamical  system: 

/*,«)][  o  Mr  |U<0  J  [  0  1 

<M‘)  J  [  —l/r  — 2/t  J  \  Ml)  J  Msfi  J  (2) 

«••,(*)  =  0i(O  +  v/50a(t), 

where  n(t)  denotes  a  stationary  white  noise  signal  that  excites  the  gust  model,  while 
and  u>j(t)  denote  the  gust  states.  We  assume  that  n(t)  has  zero  mean  and  unit  r.m.s.  value. 

That  is,  for  all  times  t,  we  have  £{n(t)}  =  0  and  E{n*(t)}  =  1.  Here,  £{•}  denotes  the 

expectation  operator.  The  time  constant  r  determines  the  dynamics  of  the  gust  model  and 
is  given  by  r  =  L/U ,  where  L  is  the  scale  length  and  U  the  airspeed.  The  model  in  Eqn. 

2  is  scaled  so  that  the  steady-state  r.m.s.  value  of  w,(t)  is  constant  and  independent  of  L. 
This  means  that  £<»{u>*(t)}  =  E{n2{t))  =  1,  where  £ao{u>2(0}  =  lim«-.oo  £{u?J(t)}  is  the 
steady-state  variance  of  u>f(t)- 

The  three  control  parameters  k,  g  and  0  are  variables  representing  the  actuator  signals 
available  for  controlling  the  structure.  The  control  input  is  given  by  the  vector 

f  MO 
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The  controls  respond  to  measured  performance  of  the  airfoil.  Let  y  denote  a  vector 
of  measurements,  corresponding  to  all  positions  and  velocities  available  for  feedback.  This 
vector  is: 

MO  ' 

y(0  =  (  1  +  MO,  (4) 

MO 

,  MO  . 

where  v(t)  represents  the  measurement  noise  vector.  We  assume  that  the  measurement  noise 
vector  is  stationary  white  noise  with  zero  mean  and  covariance  £{v(t)v(t)'}  =  VI ,  where  / 
denotes  the  identity  matrix.  The  scalar  V  is  the  noise  intensity.  It  is  also  assumed  that  the 
measurement  noise  and  the  gust  excitation  are  not  correlated.  It  should  be  noted  that  this 
information  pattern  is  not  state- feedback  in  the  usual  sense  because  the  gust  states  are  not 
directly  available  to  the  controller. 

The  airfoil  dynamics  and  the  gust  model  are  described  by  a  simple  block  diagram  shown 
in  Figure  3.  In  this  diagram  the  plant  Q  is  the  combination  of  the  structure  to  be  controlled 
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and  the  gust  model,  while  C  denotes  an  output  feedback  controller.  The  controlled  signals  are 
composed  of  performance  variables  and  scaled  actuator  variables.  The  performance  variables 
are  the  plunge  and  the  pitch  displacements,  and  they  are  defined  as: 
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The  weighting  factors  ha,  a0,  are  either  the  nominal  values  or  the  maximum  values  of  the 
corresponding  variable.  The  purpose  of  these  weighting  factors  is  to  scale  the  variables  of 
interest  so  that  they  can  be  compared  sensibly.  The  scaled  actuator  variables  are 


*•»(*)  = 


MO 

M  ’ 


*.a  (0  = 


MO 

9o 


*•3(0  = 


m 

00 


(6) 


where  the  scales,  or  weighting  factors,  k*,  g0  and  0O  are  the  maximum  values  or  the  nominal 
values  of  the  corresponding  actuator  variable. 

Active  control  of  this  model  has  two  objectives,  gust  alleviation  and  increasing  the  flutter 
speed.  For  the  gust  alleviation  problem  we  want  to  keep  the  performance  signals  z,i(0  and 
*pa(0  (MO  and  a(0)  "small”  despite  the  presence  of  the  exogenous  input  tr(0  (i-«-  the  white 
noise  n(0  driving  the  gust  model  and  the  measurement  noise  v(0)  while  also  keeping  the 
scaled  actuator  variables  *.i(0»  *.a(0.  *nd  **3(0  (M0»  MO  MO)  within  reasonable 
bounds.  We  also  want  to  use  the  active  control  to  increase  the  flutter  speed  by  at  least  10%. 

To  quantify  the  above  design  goals,  we  need  to  define  mathematically  what  we  mean 
by  a  “small”  signal.  Let  C  represent  a  linear  time  invariant  controller  that  stabilizes  the 
airfoil.  When  the  gust  defined  in  Eqn.  2  enters  the  closed  loop  system,  and  given  any  initial 
condition,  all  the  variables  in  the  system  approach  stationary  stochastic  processes  with  zero 
mean  and  r.m.s.  values  that  depend  on  the  closed  loop  system  (i.e.  they  depend  on  the  airfoil 
structure  and  the  controller).  Thus,  the  ability  and  effectiveness  of  the  closed  loop  system  to 
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reject  the  stochastic  gust  disturbance  and  the  sensor  noise  may  be  determined  by  measuring 
the  steady  state  r.m.s.  values  of  the  weighted  plunge  and  weighted  pitch  displacements. 
Similarly,  the  activity  of  each  actuator  signal  can  be  measured  by  the  steady  state  r.m.s. 
value  of  the  corresponding  weighted  signal. 

Mathematically,  given  a  plant  Q  and  a  controller  C,  the  performance  of  the  controlled 
structure  is  represented  by  the  quantities  £,*, {*£,(<)},  while  the  actuator  activity  is  given  by 
Good  disturbance  rejection  (or  good  performance)  means  small  r.m.s.  values  for 
zpl(f)  and  Zpj(t).  Thus,  the  worst  case  performance  measure  for  the  displacement  variables 
is  defined  by  the  function: 

—  m4X(^oo{*Ji(0}' ^oo{*Jj(f)})-  (7) 

The  function  takes  on  the  larger  of  two  values  of  airfoil  r.m.s.  response,  one  value 
corresponds  to  plunge  response  and  the  other  to  pitch  response. 

A  large  value  of  4>^  is  to  be  regarded  aa  “bad  performance”  because  at  least  one  of  the 
two  displacements  has  a  large  r.m.s.  value.  On  the  other  hand,  a  low  value  of  ^  is  “good 
performance”  since  both  displacements  are  jointly  small  aa  measured  by  their  r.m.s.  values. 
The  worst  case  actuator  activity  is  defined  as: 

<frvn{G,C)  =  max(£oo{z*i(f)},  £qo{*?3(0}>  £<»{*l3(0})»  (8) 

The  function  is  similar  to  dwp-  A  small  value  of  represents  a  low  cost  actuator  set 
while  a  large  value  represents  a  “high  cost”  actuator  arrangement.  It  is  reasonable  to  assume 
that  allowing  large  values  of  will  produce  small  values  of 

The  <t>  functions  introduced  in  Eqns.  7  and  8  are  called  “worst  case  figures  of  merit.”  These 
figures  of  merit  are  more  natural  and  realistic  design  parameters  than  the  usual  “weighted 
sum”  objective  function  given  by 

*,(G,C)  =  A,£oo{*,V*)}  +  (9) 

where  At  and  A3  are  positive  scalar  used  to  assign  relative  importance  to  the  individual 
objectives.  For  example,  imagine  that  the  problem  is  to  find  a  controller  C  such  that  both 
£oo{*,V<)}  and  £»{*pj(*)}  are  below  a  desirable  performance  level.  Clearly,  such  a  feasibil¬ 
ity  problem  is  mathematically  equivalent  to  the  problem  of  minimizing  (over  all  stabilizing 
controllers)  the  worst  case  figure  of  merit  4wp(0,C)  and  then  checking  to  see  if  the  minimal 
worst  case  performance  satisfies  the  given  constraint. 
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On  the  other  hand,  any  attempt  to  solve  the  feasibility  problem  by  minimizing  the  weight¬ 
ed  sum  performance  objective  vp(QX)  first  requires  the  selection  of  suitable  values  for  the 
artificial  parameters  Aj  and  A3,  which  have  little  physical  meaning.  In  addition,  the  problem 
of  weight  selection  becomes  more  complicated  if  not  impossible  when  we  add  more  perfor¬ 
mance  variables  and  we  include  the  actuator  variables  in  the  definition  of  the  weighted  sum 
performance  objective  0P(£,C). 


Tradeoff  surfaces 

To  assess  the  advantages  or  limitations  of  different  structural  parameters  (different  plants 
Q  with  different  physical  features)  for  the  problem  of  gust  alleviation  we  need  to  introduce 
a  quantitative  optimality  concept  for  the  cost  functions  and  defined  in  Eqns.  7  and 
8,  respectively. 

Since  smaller  values  of  the  figure  of  merit  functions  ^«*(£,C)  and  4>**(G,C)  are  better, 
a  plant-controller  pair  is  a  “good"  pair  if  both  worst  case  objectives  are  jointly  small.  Now, 
suppose  that  the  plant  Q  is  fixed  (the  design  of  the  airfoil  with  its  actuators  is  fixed,  but  the 
control  law  relating  actuators  output  to  sensor  input  is  not).  Then,  given  two  numbers  71 
and  7]  (representing  performance  and  actuator  r.m.s.  limits,  respectively)  we  say  that  the 
design  specification  (71,73)  is  achievable  if  we  can  find  a  stabilizing  controller  C  such  that 
we  have 

tap(G.C)  <  7i  and  4*m(G,C)  <  7a-  (10) 

From  this  equation  it  is  clear  that  the  set  of  achievable  specifications  (71,73),  which  will  be 
denoted  by  A(ff),  is  given  by  the  equation 

A(G)  =  {(71, 7a) :  there  exists  a  stabilizing  controller  C  such  that  C)  <  7i 

and  4>„{Q,C)  <  7a}-  (11) 

Although  cast  in  terms  of  control  terminology  Eqn.  11  is  central  to  the  appreciation 
of  the  power  of  this  approach.  Note  that  there  exists  at  least  one  controller  C  such  that 
the  controlled  structure  satisfies  a  given  design  specification  if  and  only  if  this  specification 
belongs  to  the  set  of  numbers  A{Q).  Note  also  that  this  set  is  a  function  of  the  airfoil/actuator 
plant  Q  only.  (This  explains  our  notation  in  Eqn.  11.)  For  instance,  if  some  of  the  structural 
parameters  are  changed,  i.e.  the  plant  Q  is  changed,  then  a  design  specification  that  was 
originally  not  achievable  could  become  achievable  with  some  stabilizing  controller. 
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MG)  defines  a  design  space  and  the  lower  boundary  of  A{G)  is  called  the  •‘tradeoff  sur¬ 
face"  or  the  set  of  Pareto  optimal  specifications  corresponding  to  the  cost  functions  0^  and 
0W«.  The  tradeoff  surface  is  a  very  important  feature  in  multiobjective  optimization  because 
it  determines  the  boundary  between  achievable  and  not  achievable  design  specifications-i.e. 
the  limits  of  achievable  performance. 

An  example  of  the  trade-off  surface  concept  to  illustrate  its  use  is  shown  in  Figure  4. 
This  figure  plots  limiting  combinations  of  71  and  7j  computed  for  two  slightly  different 
plant-actuator-sensor  combinations.  Two  curves  are  shown  in  Figure  4,  because  there  are 
two  different  plants  (controlled  structure  configurations)  to  be  compared. 


Figure  4:  Typical  trade-off  surface  (curves),  showing  feasible  and  nonfeasible  regions 
for  two  different  plants. 

For  each  plant,  combinations  of  worst  case  performance  (71)  and  worst  case  actuator 
cost  (7a)  lying  below  each  curved  line  are  not  possible  unless  some  feature  of  the  system  is 
changed.  This  change  might  be  relocation  of  an  actuator,  resizing  of  an  actuator,  changing 
the  stiffness  of  the  airfoil  or  changing  the  actuator  material.  Note  that  these  two  trade-off 
curves  intersect.  This  means  that,  for  one  range  of  performance  requirements,  one  plant 
configuration  is  better  than  the  other. 

Trade-off  surfaces  are  also  valuable  once  control  laws  are  proposed.  The  tradeoff  surface 
may  be  used  as  an  absolute  measure  against  which  different  control  laws  can  be  compared. 
Thus,  the  limits  of  achievable  performance  are  valuable  not  only  to  determine  which  struc- 
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tural  configuration  is  easier  to  control  but  also  they  provide  the  necessary  information  to 
evaluate  candidate  controller  designs  with  respect  to  the  figures  of  merit  of  interest. 

The  methodology  to  compute  trade-off  surfaces  (in  our  case  these  are  2-dimensional  so 
they  become  trade-off  curves)  for  the  figures  of  merit  ^  and  <t>^  is  taken  from  [4]  and  [5], 
For  an  e:  v  >sition  of  the  area  of  multiobjective  control  see  (3).  See  also  [6]  for  some  more 
recent  work  in  this  field.  It  is  important  to  note  that  these  trade-off  surfaces  are  computed 
under  the  only  assumption  that  the  controllers  to  be  used  are  linear  and  time- invariant.  No 
other  restriction  on  the  structure  of  the  controllers  (i.e.  the  order  of  the  controllers  should 
be  bounded,  the  controllers  should  have  a  specific  architecture  such  as  observer-based,  etc.) 
are  imposed.  This  means  that  if  a  design  specification  (71,73)  lies  below  a  curved  line  in 
Figure  3,  there  is  no  stabilizing  controller  (even  infinite-dimensional)  that  can  achieve  this 
specification. 

While  the  mathematical  theory  to  create  trade-off  surfaces  haa  been  well  developed  in 
the  above  references,  the  programming  of  the  numerical  procedure  requires  time  and  skill.  A 
preliminary  version  of  the  software  utilized  in  this  work  may  be  found  in  [5].  This  software 
makes  use  of  commercial  tools  for  control  systems  analysis  and  design,  and  finite  dimensional 
convex  optimization  algorithms.  This  preliminary  software  package  has  been  implemented 
in  Matlab.  The  numerical  methods  used  not  only  provide,  to  any  desired  degree  of  accuracy, 
points  in  the  trade-off  surface  but  also  controllers  whose  cost  values  are  as  close  as  desired 
to  the  trade-off  surface. 


Numerical  examples 

In  this  section  a  tradeoff  surface  between  the  performance  cost  (measured  in  terms  of 
worst  case  r.m.s.  airfoil  response  in  pitch  or  plunge)  and  the  actuator  cost  (measured  in 
terms  of  its  r.m.s.  value)  is  computed  for  three  different  cases.  In  each  case,  the  plants 
use  a  different  actuator.  This  simple  problem  will  illustrate  the  features  and  value  of  our 
approach.  We  will  conclude  that  the  torsional  spring  actuator  g  is  the  best  actuator  signal, 
if  only  one  actuator  can  be  used,  for  the  gust  alleviation  problem  at  an  airspeed  larger  than 
the  open-loop  flutter  speed. 

The  closed-loop  design  airspeed  is  set  to  10%  above  the  airfoil  uncontrolled  flutter  speed. 
The  weights  for  the  performance  variables  z,  and  the  actuator  variables  za  defined  in  Eqns. 


16 


1 


differently.  For  the  bending  strain  actuator  k  (daahed)  we  plot  vs.  =  ^{^(t)} 
For  the  torsion  strain  control  g  (solid)  we  plot  vs.  10  4**  =  10  £«>{^i(t)}.  Finally,  for 
the  trailing  edge  flap  control  (dashdot)  we  plot  vs.  0.1^*,  =  0.1  £oo{*2j(0}-  The  torsion 
strain  actuator  g  is  seen  to  be  the  moat  effective  for  all  possible  performance  values. 

The  second  example  uses  a  gust  scale  length  L  =  10  b  but  keeps  all  other  parameters 
fixed.  This  results  in  the  Dryden  transfer  function  having  bandwidth  on  the  order  of  the 
open  loop  torsion  resonant  frequency.  For  example,  r  =  0.07726  sec.  and  u>*r  =  3.863.  so 
that  uja  is  of  the  same  order  of  the  corner  frequency  of  the  Dryden  transfer  function. 

Figure  6  shows  the  tradeoff  surface  between  the  worst  case  performance  index  ^  and 
the  worst  case  actuator  index  for  the  torsion  and  bending  strain  controls  g  and  k.  As 
before,  we  plot  ^  vs.  =  £<x,{*2i(0}  for  the  bending  strain  actuator  k  (dashed)  and 
4)^  vs.  10 =  10£,»{*ia(O}  for  the  torsion  strain  control  g  (solid).  Clearly,  torsion  strain 
control  is  more  effective  than  bending  strain  control. 
sio4 


Figure  6:  Trade-off  surface  for  active  material  bending  and  torsion  control. 

Gust  scale  =  10  b. 

Finally,  in  Figure  7,  we  compare  trailing  edge  flap  control  (dashdot)  with  bending  strain 
control  (dashed).  This  plot  does  not  use  scaling  so  that  actual  costs  are  shown.  It  is  again 
obvious  that  the  active  strain  control  is  better  than  the  trailing  edge  flap  control  in  these 
cases. 
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Aileron  control  Actuator 


Bonding  strain  actuator 


"o  o.ooi  o.oo  0.003  um  tom  <xm  amt  un  ooo*  oni 

Worst  cast  actuator  cost 

Figure  7:  Trade-off  surface  for  active  beading  control  and  aileron  control. 

Future  efforts  and  summary 

The  work  thus  far  has  concentrated  on  developing  a  pilot  program  method  to  evaluate 
actuator /structure  aeroelastic  control  effectiveness  so  that  the  effectiveness  of  integrated 
active  ma'  -nal  actuators  can  be  predicted.  By  any  measure,  this  effort  has  been  a  success. 
However,  if  .ae  method  is  to  have  an  impact  of  the  controlled  structure  community,  it  will 
have  to  be  demonstrated  for  a  more  elaborate  analytical  model  where  additional  design 
parameters  can  be  exposed. 

At  the  present  time  a  new  aeroelastic  model  is  being  developed  to  simulate  a  closed- 
cell  box  beam  structure  with  embedded  piezoelectric  elements.  This  finite  element  model 
contains  unsteady  aerodynamic  lags  and  actuator  dynamics.  The  finite  element  model  will  be 
more  useful  in  the  assessment  of  actual  active  material  capabilities  and  actuator  placement 
and  sizing.  Additions  to  the  present  trade-off  surface  software  and  the  study  of  trade-off 
surfaces  for  other  figures  of  merit  of  practical  significance  are  necessary.  Other  figures  of 
merit  we  plan  to  consider  include  the  maximum  time-domain  excursion  of  the  signals  of 
interest  in  response  to  fixed  or  worst  case  exogenous  inputs. 

The  results  of  this  work  will  be  used  to  study  active  material  capabilities,  assess  directions 
for  improvement  in  this  area,  study  sizing  and  placement  are  create  data  that  can  be  used 
for  design  of  future  test  articles. 
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ABSTRACT 

This  papw  discusses  the  development  and  use  of  an  effective  finite  element  analysis  procedure  id  examine  integrated 
anisotropic  piezoelectric  actuators  for  aeroelastic  applications  where  it  is  essential  that  torsion  and  bending  be  controlled 
independently  of  one  another.  The  finite  element  model  consists  of  three-dimensional  isoparametric  solid  elements  (bricks) 
that  allow  modeling  of  tailored  piezoelectrics  with  skewed  actuator/sensor  axes.  This  scheme  also  allows  the  representation 
of  an  anisotropic  host  structure  and  can  account  for  material  and  stacking  geometry  through  the  element  thickness.  Using  this 
finite  element  representation,  it  is  shown  that  anisotropic  piezoelectric  actuators  can  create  sufficient  twisting  and  bending  to 
control  aerodynamic  loads  on  a  wing,  although  aerodynamic  loads  are  not  included  in  this  discussion. 

1.0  BACKGROUND  AND  PURPOSE 

The  purpose  of  this  paper  is  to  examine  the  action  and  effectiveness  of  adaptive  material  actuators  when  they  ate 
used  to  create  deflections  of  wing  structures.  By  effectiveness,  we  men  the  ability  id  cause  twist  and  bending  deflections 
that  will,  in  turn,  create  aerodynamic  loads.  Here  the  emphasis  will  be  on  structural  deflection,  not  aerodynamic  interactions, 
so  no  aerodynamic  loads  will  be  included. 

Adaptive  structures  use  actuators  to  create  changes  in  structural  design  stupe  or  damping  to  respond  to  dunging 
mission  needs  and  performance  requirements.  Actuators  may  be  articulated  mechanical  devices  such  as  ailerons  and  flaps  or 
they  may  be  'smart''  materials  embedded  in  the  structure.  In  the  case  of  atmospheric  flight  vehicles,  active  material  actuators 
embedded  in  the  structure  can  control  structural  shape  and  aerodynamic  loads.  Using  sensors  and  feedback  control 
algorithms,  these  actuators  also  can  change  wing  stiffness  and  control  the  stability  of  an  entire  lifting  surface  to  maneuver  and 
trim  an  aircraft,  reduce  gust  loads  or  enhance  the  stability  of  the  structure.  In  some  cases,  this  can  be  done  with  less  weight 
and  at  reduced  overall  cost  compared  to  conventional  structures.  At  present,  these  concepts  rely  on  aeroeiasticity.  that  is,  the 
mutual  interaction  between  aerodynamic  loads  and  structural  deflections,  to  create  favorable  active  load  control, 

Aeroelastic  control  concepts  are  still  in  their  infancy,  but  their  proposed  uses  are  increasing  rapidly  with  each 
passing  year.  The  successful  application  of  active  materials  and  actuator  concepts  depends  on:  a)  identification  of  an 
aeroelastic  phenomenon  and.  (b)  identification  of  an  effective  actuator  arrangement  to  control  the  phenomenon.  Both  involve 
creativity  and  the  accessibility  of  an  effective  analytical  procedure.  The  use  of  actuators  to  control  aeroelastic  deformation 
also  depends  on  their  ability  to  decouple  two  characteristic  types  of  deformation,  bending  and  twist  This  modal  decoupling 
is  very  important  to  aeroelastic  load  control  because  both  types  of  deformation  create  aerodynamic  loads,  but  we  may  not 
always  want  these  loads  to  occur  simultaneously. 

The  study  of  aeroelastic  decoupling  requires  an  accurate  analytical  model,  one  that  can  model  arrays  of  anisotropic 
actuators  oriented  at  angles  skewed  with  respect  to  conventional  structural  axes.  It  also  requires  realistic  modeling  of  an 
aircraft  wing  structure.  This  model  may  be  a  plate-like  configuration  or  the  closed  box-beam  arrangement  common  to  all 
efficient  aeronautical  structures. 

Because  of  the  interest  in  effective  actuator  combinations  and  realistic  structures,  a  finite  element  method  was 
developed  to  model  potential  structure/actuator  configurations.  This  technique  draws  heavily  upon  previous  published  work 
by  several  researchers,  but  it  extends  this  work  to  cover  emerging  structural  configurations.  This  theoretical  development 
will  be  reviewed  and  discussed  in  this  paper. 
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We  first  turn  our  attention  to  a  survey  of  previous  active  structures  work  in  the  area  of  aeroeLasoc  and  aerodynamic 
control.  This  review  will  help  the  reader  to  understand  the  types  of  problems  that  are  potentially  worthwhile  to  pursue  and 
the  difficulties  applying  these  concepts.  Having  done  this  we  summarize  the  finite  element  development  and  discuss 
illustrative  examples. 

1.1  Aeroelastlc  and  aerodynamic  load  control 

Control  of  aerodynamic  loads  and  structural  response  of  wings  and  other  lifting  surfaces  is  not  a  new  concept. 
Proposals  to  control  wing  response,  in  particular  gust  response  are  over  thirty  years  old.  Ride  quality  enhancement  was  used 
on  the  B-70  VaUcerie  Supersonic  Bomber.  The  same  concept  is  used  today  on  the  B-l  bomber,  to  allow  it  to  conduct  low 
level .  high  speed  penetration  missions.  The  Lockheed  L-101 1  uses  active  inboard  flaps  to  reduce  the  severity  of  turbulence 
on  wing  fatigue  life  and  on  the  ride  quality  felt  by  passengers.  Suppression  of  wing  flutter,  a  dynamic  structural  instability 
created  by  unfavorable  interaction  between  unsteady  aerodynamic  loads  and  structural  vibrations,  was  demonstrated  in  the 
mid- 1970's.  This  flutter  suppression  was  possible  by  using  feedback  control  to  drive  aerodynamic  surfaces  such  as  ailerons 
and  leading  edge  slats  to  damp  out  motion. 

Active  materials,  in  particular  piezoelectric  material  actuators,  have  been  proposed  as  controllers  to  accomplish 
much  the  same  mission  as  aerodynamic  surfaces.  How  successful  and  in  what  areas  these  material  actuators  will  be  able  to 
replace  or  augment  aerodynamic  controllers  such  as  ailerons  is  still  very  much  in  doubt  There  are  advantages  to  creating  so- 
called  "solid  state  wings"  that  have  no  external  hinges  or  mechanical  devices.  For  military  applications  there  is  the  advantage 
of  stealth.  Any  gaps  or  openings  that  appear  in  the  wing  tends  to  increase  size  of  the  radar  return  from  an  aircraft  Another 
advantage  is  that  the  control  is  distributed  over  the  entire  wing,  making  control  of  some  aeroelastic  phenomena  more  precise. 
Lastly,  there  are  some  aeroelastic  phenomena,  such  as  panel  flutter,  for  which  no  mechanical  device  is  feasible. 

The  first  attempt  at  aerodynamic  shape  control  is  due  to  Crawley,  etal\  at  M.I.T.  They  examined  the  effectiveness 
of  piezoelectric  actuators  to  create  twist  and  camber  in  plate-like  aerodynamic  surfaces.  The  ability  of  these  actuators  to 
affect  aeroelastic  loads  was  analyzed  and  demonstrated  in  wind  tunnel  experiments  by  Lazarus  and  Crawley2.  In  a  related 
study,  Bohlmann  and  Lazarus3  examined  the  response  of  piezoelec Uically  actuated  flat  plates  in  an  airstream.  For  this  latter 
study,  piezoelectric  actuation  created  chord  wise  plate  curvature  (which  creates  an  equivalent  angle  of  attack)  to  generate 
aerodynamic  forces. 

While  early  studies  emphasized  the  ability  of  piezoelectric  actuators  only  to  create  airloads,  the  essential  aeroelastic 
feedback  process  was  not  well  understood.  At  Purdue  University,  Ehlers  and  Weiss  ha*4-5’6  formulated  a  simple  aeroelastic 
model  for  a  beam-like  laminated  composite  wing  with  embedded  actuators.  This  model  captured  the  fundamental 
interactions  between  bending  and  torsional  deformations  and  aerodynamic  loads.  Using  this  model,  they  showed  that  an 
aeroelastic  static  instability  called  divergence  could  be  controlled.  They  also  demonstrated  the  wing  lift  effectiveness  (the 
ability  to  increase  or  decrease  lift)  could  be  changed  with  a  simple  feedback  control  law  relating  wing  root  bending  and 
twisting  moments  and  piezoelectric  action.  They  developed  aero-piezo-elasdc  parameters  for  measuring  the  effectiveness  of 
actuators  in  creating  aerodynamic  loads. 

Lazarus,  et  ai.  7  have  shown  the  potential  benefits  that  may  occur  when  aerodynamic  control  surfaces  are  replaced 
with  strain  actuation.  These  results  are  limited  to  simple  plate  configurations.  Ehlers4  has  also  studied  ways  to  increase 
aileron  effectiveness  using  a  combination  of  strain  actuation  and  ailerons. 

Song,  et  at.  ®  considered  static  aeroelastic  control  using  a  thin-walled  beam  structure  to  represent  a  wing  structure. 
This  formulation  is  important  to  aeroelastic  control  analysis  because  the  box-beam  is  the  preferred  configuration  for  wing 
structural  design,  although  plate  models  may  capture  most  of  the  wing  behavior  when  the  wing  is  thin.  Persiani,  et  ai.  ?  have 
al«n  considered  the  static  aeroelastic  behavior  of  an  adaptive  structural  box  beam  with  composite  materials.  Their  study  used 
the  concept  of  strain  energy  tuning. 

Dynamic  aeroelastic  phenomena  have  also  received  attention.  Lazarus,  Crawley  and  Lin*®  showed  that 
piezoelectric  strain  actuation  might  be  a  viable  alternative  to  articulated  control  surfaces  for  flutter  suppression.  Panel  flutter 
suppression  has  been  examined  by  Scott*  *,  Abou-Amer*2,  Paige*3,  Hajela  and  Glowasky*4  and  Scott  and  Weisshaar*5. 
Panel  flutter  suppression  with  piezoelectric  actuators  and  shape  memory  alloy  actuators  is  unique  in  that  no  articulated  device 
exists  to  do  the  same  task.  Heeg  *^  has  also  demonstrated  airfoil  flutter  control  with  a  piezoelectric  actuator  device. 
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In  all  previous  results,  shape  control  of  aerodynamic  surfaces  was  based  on  the  inherent  extension/contnction  ability 
of  piezoelectric  (usually  piezoceramic)  actuator  materials.  This  action  was  used  to  create  bending  distortion  of  plate-like  or 
beam -like  materials,  which,  in  turn,  tend  to  twist  the  structure.  This  twist  was  made  possible  only  if  there  was 
bending/torsion  coupling  in  the  wing  structure,  such  as  will  occur  if  laminated  composites  are  used  or  if  there  is  geometrical 
bend/twist  coupling  in  swept  wing  planforms.17-18  While  some  of  these  previous  studies  accounted  for  the  possibility  of 
creating  shear  to  cause  torsion,  it  is  generally  recognized  that  the  most  popular  piezoceramics  are  electromechanically 
isotropic  so  that  shear  strain  generation  is  not  possible. 

Obtaining  decoupled,  direct  control  of  wing  twist  is  highly  desired  for  aeroelasdc  and  aerodynamic  load  generation. 
Innovative,  alternative  ways  to  create  decoupled  twist  have  been  suggested  by  Lee19  and  Barrett.20-21  Barrett's  proposed 
configuration  uses  a  directional  attachment  of  piezoelectric  strips  (DAP)  to  channel  the  actuation  strain  into  torsion. 
Reference  20  is  among  a  rapidly  increasing  number  of  applications  to  rotorcraft  noise  and  vibration  reduction.  Although  this 
procedure  has  been  demonstrated  to  be  effective  in  wing  tunnel  experiments,  the  complexity  of  the  DAP  construction  may 
reduce  its  effectiveness  for  actual  applications. 

12  Objectives  -  Independent  bending/tonioa  control  and  finite  element  modeling 

One  objective  of  this  paper  is  to  explore  alternative  ways  of  obtaining  direct  actuation  control  of  adaptive  structures 
to  create  a  desired  mode  of  wing  twist  and  bending.  An  active  tailoring  concept  based  on  active  orthotropicity  of  some  types 
of  piezopolymeric  materials  will  be  introduced.  Until  now,  these  piezopolymerics  have  been  excluded  from  aeroelasdc 
control  applications  because  they  are  not  stiff  and  because  they  require  high  voltages  to  create  the  electric  fields  necessary  to 
be  effective.  However,  if  the  applied  voltage  is  not  limited,  an  acceptable  value  of  the  piezoelectric  effectiveness  parameter 
defined  by  Crawley  and  DeLuis22  can  be  obtained. 

Die  unique  property  of  electromechanical  anisotropy  in  special  kinds  of  PVDF  piezopolymers  is  obtained  by 
stretching  during  poling.23*24,  This  property  was  used  for  the  first  time  by  Lee  and  Moon25  to  design  micro  actuators  and 
sensors.  However,  applications  of  PVDF  to  adaptive  structures  have  been  limited  to  shape  control  and  its  use  in  structural 
damping2**-27  and  in  pointing  devices.28-29-30  In  all  of  these  applications,  "ordinary"  PVDF  with  isotropic 
electromechanical  properties  was  used. 

Decoupled  piezoelectric  control  and  the  ideas  related  to  such  control  require  an  effective  analytical  model. 
Developing  such  a  model  is  the  second  objective  of  this  study.  Structural  models  of  structures  with  integrated  actuators  and 
sensors,  bonded  or  embedded  in  the  host  structure,  have  bean  limited  mostly  to  simple  beam  models  or  rectangular  plates. 
Beam  models  have  permitted  closed-form  solutions  to  some  important,  problems.22-31 

The  Rayleigh-Ritz  method  has  also  been  applied  to  active  structures.32  This  method,  when  used  in  conjunction  with 
plate  theory  or  beam  theory  can  be  effective.  On  tiie  other  hand,  this  procedure  requires  a  careful  choice  of  approximating 
functions  and  is  difficult  to  apply  to  general  cases  with  arbitrarily  shaped  actuators  and  structural  planfocms. 

Finite  element  models  have  also  been  developed  to  analyze  active  structures.  These  efforts  have  fallen  into  two 
categories:  (a)  modeling  pure  piezoelectric  structures  for  transducer/vibration  studies33-34*35-36-37;  and  (b)  modeling  host 
structures  with  integrated  piezoelectric  actuators/sensors,  bonded  or  embedded  in  their  host  structures.  38 -39-40  For  these 
latter  studies,  the  finite  element  method  was  used  only  on  rectangular  plates,  not  closed  thin-wall  sections.  More  importantly, 
these  studies  also  were  restricted  to  piezoelectric  actuators  with  isotropic  electromechanical  properties. 

This  paper  outlines  the  development  of  a  finite  element  model  developed  to  analyze  plate  based  structures  (including 
box-beams)  subjected  to  mechanical/electrical  excitations.  The  finite  element  model  uses  a  thin  piezoelectric  element  derived 
from  Hamilton's  Principle.  The  theory  is  restricted  to  linear  piezoelectric  theory,  but  takes  into  account  orthotropic 
mechanical  and  electromechanical  properties  of  the  actuators. 

2.0  FINITE  ELEMENT  DEVELOPMENT 

The  purpose  of  this  section  is  to  review  the  development  of  the  finite  element  model  used  to  analyze  integrated, 
active,  piezoelec tri call y  activated  structures.  This  review  will  include  development  and  restrictions  on  its  use. 
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A  piezoelectric  material  is  a  material  which,  when  subjected  to  mechanical  load,  develops  an  e lectnc  charge 
proportional  to  the  resulting  mechanical  stress.  Conversely,  this  material  will  deform  or  strain  when  an  electric  Held  is 
applied.  This  deformation,  in  particular  positive  or  negative  strain,  depends  on  the  polarity  of  the  applied  field.  The  electric 
field  is  generated  by  imposing  a  voltage  across  the  material  so  that  the  field  is  proportional  to  the  applied  voltage  divided  by 
the  distance  between  the  electrodes. 

Figure  1  (all  figures  have  been  placed  at  the  end  of  this  paper)  represents  the  general  type  of  structure  to  be  modeled. 
The  structure  is  composed  of  an  arbitrary  plate  with  embedded  discrete  piezoelectric  elements.  These  actuator  elements  are 
poled,  elec  traded  and  connected  arbitrarily.  The  equations  of  motion  of  this  continuum  can  be  derived  using  a  generalized 
form  of  Hamilton's  principle.  The  variational  expression  for  this  coupled  electromechanical  system  can  be  written  as  follows: 

tf,2[5(r  ~U~  Wdtc  )  +  Swea  ]dr  =  0  (1) 


In  Eqn.  1  the  following  definitions  are  used:  T  =  kinetic  energy ;  U  =  poundal(strain)  energy, 

Welgc  =  electrical  energy ;  SWgja  =  virtual  work  of  external  forces . 

While  Jk  expressions  for  the  structural  kinetic  energy  and  strain  energy  are  the  standard  expressions  for  a  continuum,  the 
electrical  energy  and  the  virtual  work  from  the  external  forces  include  piezoelectric  constitutive  relations  and  properties. 

A  piezoelectric  element  is  shown  in  Figure  2,  together  with  its  electrodes  and  coordinate  axis  definitions.  The 
general  properties  of  a  piezoelectric  material  are  defined  relative  to  the  material  poling  direction.  This  poling  direction  is 
taken  parallel  to  the  X3  direction  .  We  require  a  linear  relationship  between  the  elastic,  electrical  and  electromechanical 

coupling  properties  of  the  piezoelectric  material.  This  relationship  is  given  by  its  piezoelectric-dielectric  relations.41*42 
These  relations  read  as  follows: 


(2) 


Equation  2  relates  the  electrical  displacement  (charge)  D  and  the  mechanical  stress  a  to  the  applied  electric  field,  E,  and  the 
mechanical  strain  S  through  constants  of  proportionality  representing:  dielectric  constants  e  for  the  material;  e,  piezoelectric 

constants  relating  voltage  to  stress;  and  Cp,  the  piezoelectric  material  stiffness  matrix  elements.  The  notation  (  ) 

/ 

indicates  that  the  constants  are  measured  at  constant  strain  (clamped),  while  (  )  indicates  values  computed  in  the  local 

piezoelectric  material  coordinate  reference  frame.  The  superscript  notation  (  )£  indicates  that  values  are  measured  at 
constant  electric  field  (short  circuited). 

The  relationship  in  Eqn.  2  can  be  inverted,  either  in  its  entirety  or  partially.  In  particular,  the  following  relationship 
is  valuable. 

KI-kT'M  ® 


The  elements  d-  are  the  piezoelectric  charge  coefficients  that  relate  induced  strain  to  the  applied  electric  field.  Equation  3 
can  be  used  to  calculate  the  etj  constants  in  terms  of  the  d^  constants,  since  the  latter  constants  are  generally  more  available 
than  the  constants. 

2.1  Typical  piezoelectric  materials  and  properties 


Table  1  summarizes  some  important  elastic  and  electromechanical  constants  for  lead  zirconate  titanate  (PZT)  and 
polyvinylidene  fluoride  (PVDF).  An  examination  of  the  parameter  values  in  this  table  reveals  some  important,  essential 
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features  of  these  materials.  Fust  of  all,  note  the  signs  of  the  parameters  d,,  in  Table  1 .  If  these  coefficients  are  negative,  the 
voltage  must  be  applied  in  the  direction  opposite  to  the  poling  direction  to  produce  positive  (extensional)  strains.  Note  also 
that  the  induced  strain  in  PVDF  is  nearly  uniaxial. 

Piezoceramics  have  two  orthogonal  planes  of  symmetry.,  one  inpiane  and  the  other  out-of-plane  so  that,  while  they 
are  orthotropic  materials,  there  is  planar  isotropy  in  the  plane  perpendicular  to  the  poling  axis.  This  planar  isotropy  is  evident 
in  Table  1  for  PZT  where  d32sd3i  and  d24=dij.  On  the  other  hand,  piezopolymencs  such  as  PVDF  have  fully  isotropic 
structural  properties,  but  electrically  induced  inplane  strains  that  are  decidedly  orthotropic.  Because  the  orthotropicity  is  in  a 
plane  perpendicular  to  the  poling  axis,  it  presents  a  potential  for  actuator  tailoring. 

Tahiti 

Piezoelectric  material  constants  for  PZT  and  PVDF 


units 

PVDF  [241  S 

El 1*E22 

GPa 

0.63xl0u 

£33 

GPa 

0.49xl0n 

<331=032 

GPa 

0-22x10** 

Gj2 

GPa 

0.233x10** 

0.77xl09 

V12 

0.35 

0.3 

d31 

m/V 

-166-xlO**2 

0.23xl0**0 

d32 

m/V 

-166-XlO**2 

0.3x10* 1 

d33 

m/V 

360.X10*12 

-0.33x10**° 

d24 

m/V 

540.X10-12 

0 

<*15 

m/V 

540.X10'12 

0 

F/M 

15-xlO9 

0.1062xl0*9 

22  Inducing  torsion  with  skewed  actuators  and  actuator  material  axes 

When  the  principal  strain  axes  of  orthotropic  piezopolymers  are  rotated,  or  skewed,  with  respect  to  the  host 
structure,  we  can  induce  shear  strains.  In  plates  or  beams,  this  shear  creates  twisting  of  the  structure.  Figure  3  illustrates  this 
behavior  on  a  simple  skewed  plate.  A  PVDF  actuator,  skewed  at  an  angle  a,  creates  a  shear  stress  given  as: 

t  =  <r12  ■  ±V3  smacosa(*3}  -  e32)  (4) 

This  shear  stress  creates  an  equivalent  shear  flow  in  the  plate,  given  as  =  «  =  a\2l-  Note  that,  if  the  d3i  and  d32 
coefficients  an  equal,  then  (63  j  •  C32)  is  zero  and  the  twisting  shear  flow  will  be  zero. 

Orthotropic  action  of  actuators  is  essential  far  any  successful  aeroelastic  control  scheme.  The  twisting  action  of  the 
actuator  comes  from  two  sources:  (a)  resultant,  effective  torque  at  the  free  end  of  the  plate  created  by  opposing  shear  flows  on 
the  upper  and  lower  actuator  elements;  and,  (b)  differential  bending  applied  to  the  two  parallel  free  edges  of  the  plate.  The 
differential  bending  induces  warping  in  the  ’correct”  direction  of  twist 

23  The  elemental  model 

The  basic  building  block  finite  element  model  used  in  this  study  is  a  solid  brick  element  shown  in  Figure  4.  The 
degrees  of  freedom  of  this  element  include  translational  displacements  and  electrical  displacements  at  each  node.  The 
electrical  displacements  correspond  to  the  electrical  potential  4.  There  are  three  mechanical  (translational)  degrees  of 
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freedom  and  one  electrical  degree  of  freedom  per  node  and  there  are  8  nodes  per  element,  as  indicated  in  Figure  4.  This 
element  has  a  trapezoidal  planform,  with  constant  thickness. 


The  shape  functions  for  the  mechanical  displacement  and  electrical  potential  Helds  are  found  in  References  44  and 
45.  The  displacement  Held  includes  degrees  of  freedom  for  incompatible  modes  using  a  'bubble  function”  found  in 
Reference  46.  These  incompatible  displacements  and  modes  were  added  to  avoid  element  locking  under  parasitic  shear  and 
to  soften  the  extra  thin  brick  element  in  the  transverse  direction.  These  shape  functions  were  used  to  generate  elemental 
mass,  stiffness  and  external  load  matrices  that  were  then  assembled  to  obtain  a  structural  equilibrium  equation  for  the 
discretized  structure.  This  equation  has  the  following  form: 


where  u  represents  mechanical  nodal  displacements  and  4  represents  electrical  displacements  related  to  sensor  output  from 

the  piezoelectric  elements.  The  notation  u  and  4  represent  differentiation  with  respect  to  time.  The  Mij  elements  are 
inertia  elements,  while  the  Kjj  elements  are  stiffness  matrix  elements.  Equation  5  shows  that  there  is  electro-mechanical 
coupling  in  the  structure.  Note  that  no  electrical  "inertia*  or  damping  is  included  in  this  model. 

This  finite  element  procedure  was  incorporated  into  a  computer  code  that  can  analyze  both  flat  plates  and  box -beams 
with  rectangular  cross-sections.  This  computer  code  and  extensive  details  on  the  theory  are  included  in  Reference  47.  This 
code  is  able  to  accurately  model  a  variety  of  wing  structural  models,  ranging  from  plates  to  beams  to  thin-wall  built  up 
structures  with  stiffeners.  In  the  examples  to  follow,  we  will  demonstrate  some  of  these  capabilities,  together  with  some 
essential  features  of  anisotropic  actuators. 

3.0  RESULTS  AND  DISCUSSION 


To  illustrate  the  application  of  the  finite  element  model  described  in  the  Section  2.0,  we  have  applied  this  method  to 
four  different  cases.  These  cases  illustrate  the  orthotropic  actuator  use  and  response.  Case  1  examines  a  bimorph  plate 
constructed  of  two  PVDF  layers  so  that  the  top  layer  acts  as  a  distributed  actuator  while  the  bottom  layer  acts  as  a  sensor. 
Case  2  is  a  rectangular  aluminum  plate  sandwiched  between  two  piezoceramic  layers.  Case  3  considers  the  same  aluminum 
plate,  but  with  two  layers  of  tailored  piezopolymers.  Finally,  Case  4  examines  a  steel  plate  sandwiched  between  two  PVDF 
elements  with  off-axis,  mirror  symmetry. 

Consider  Case  1  with  the  1  mm.  thick  bimorph  cantilever  beam  constructed  of  two  PVDF  layers  (with 
electroAnechanical  properties  given  in  Table  1),  illustrated  in  Figure  5.  The  lop  layer  acts  as  an  actuator  and  the  bottom  layer 
acts  as  a  sensor.  As  a  result,  the  sensor  signal  (nodal  electrical  potential  4s)  »*  proportional  to  the  induced  mechanical  strain 
in  the  bottom  layer.  This  beam  is  subdivided  into  20  identical  finite  elements.  The  major  strain  axis  is  aligned  in  the  xj 
direction  for  both  top  and  bottom  layers  so  dial  the  skew  angle  is  90°  so  that  no  anisotropic  coupling  is  created. 

Figure  6  shows  the  beam  deflection  and  sensor  output  due  to  an  electrical  excitation  of  1000  volts,  applied  across  the 
top  actuator  layer.  Because  the  actuator  extends  uniformly  from  the  cantilever  beam  root  to  its  tip,  the  effect  of  the  actuator 
on  the  beam  is  the  same  as  if  a  concentrated  moment  were  applied  at  the  tip.  As  a  result,  uniform  bending,  with  a  constant 
strain  level  occurs,  as  indicated  by  the  nearly  uniform  sensor  output  shown  in  Figure  6.  Only  the  beam  root  region,  where 
the  finite  element  method  correctly  predicts  a  boundary  strain  transition  region,  is  the  strain  not  uniform. 

Case  2,  considers  a  1  mm.  thick  uniform  aluminum  plate  sandwiched  between  two  0.25  mm.  piezoceramic  layers, 
with  asymmetric  poling,  as  shown  in  Figure  7.  Sixty  finite  elements  are  used  to  capture  the  plate-like  behavior  when  the  two 
actuator  layers  are  activated  by  100  Volts,  applied  asymmetrically.  As  a  result,  one  actuator  stretches  while  the  other 
contracts.  As  shown  in  Figure  8,  the  deflection  of  the  plate  is  composed  primarily  of  span  wise  bending.  However,  some 
camber  bending  occur  because  of  the  influence  of  dn-  It  is  impossible  to  generate  twist  with  this  actuator  pair  arranged  in 
this  way. 


In  an  attempt  to  generate  twist ,  the  same  host  (date  is  activated  with  offset  (skewed)  piezoceramic  elements,  as 
shown  in  Figure  9.  The  resulting  displacement,  shown  in  Figure  10.  is  that  a  considerable  amount  of  twist  is  generated,  but 
substantial  downward  bending  also  accompanies  this  twist 


26 


Case  3  uses  the  same  host  structure/actuator  arrangement  as  shown  in  Figure  7,  with  the  exception  that  the  actuators 
are  0.25  mm.  thick  PVDF.  These  onhotroptc  PVDF  layers  are  arranged  with  their  major  axes  at  45°  on  the  upper  side  and 
-45°  on  the  lower  side  to  achieve  mirror  symmetry  of  the  actuator  pair.  When  10,000  Volts  are  applied  to  each  nrtnat^r  layer, 
the  resulting  displacement,  shown  in  Figure  11,  is  pure  twist  (except  in  the  boundary  region  near  the  root).  In  addition,  the 
twisted  sections  do  not  have  any  chord  wise  deformation  in  the  chord  wise  (camber)  direction. 

Case  4  considers  a  steel  plate  with  PVDF  actuators  attached,  as  illustrated  in  Figure  12.  The  principal  coordinate 
axes  of  the  PVDF  actuators  are  arranged  at  45°  on  the  top  and  -45°  on  the  bottom.  This  configuration  corresponds  to  one 
used  by  Barren  in  Reference  21.  A  voltage  of  32,000  volts  is  then  applied.  This  voltage  is  large,  but  it  represents  the 
maximum  allowable  voltage  per  unit  thickness  (40  Volts/pm)  that  is  given  by  the  manufacturer.  The  resulting  torsional 
displacement,  plotted  as  afuncdon  of  distance  from  the  cantilever  support,  is  shown  in  Figure  13. 

The  tip  displacement  shown  in  Figure  13  is  5.25  degrees  and  is  comparable  to  results  obtained  by  Barren  for  his 
EDAP  construction.  Despite  the  very  large  voltage  used  here,  the  simplicity  of  the  PVDF  technology  may  make  it  a 
candidate  for  twister  actuators. 

It  is  important  to  note  the  influence  of  the  thickness  of  the  PVDF  layer.  Increasing  the  actuator  thickness  will 
increase  the  induced  shear  flow  and  the  induced  torque,  but  it  will  also  increase  the  torsional  stiffness  of  the  section.  As  a 
result,  there  will  be  an  optimal  thickness  for  maximum  twist  In  the  present  case,  additional  calculations  showed  this  optimal 
actuator  thickness  to  be  about  0.8  mm. 


4.0  CONCLUSION 

An  effective  and  versatile  method  of  finite  element  analysis  to  investigate  adaptive  structure  combinations  with 
actuators  having  in-plane  electromechanical  orthotropicity  has  been  outlined.  This  analytical  procedure  was  used  to 
demonstrate  the  ability  of  piezoelectric  actuators  to  generate  pure,  decoupled  torsion  for  aeroelastic  control.  Among  the 
conclusions  of  this  study  are: 

a)  Electromechanical  anisotropy  is  an  effective  and  highly  desirable,  if  not  essential,  to  aeroelastic  control. 
Efforts  should  be  directed  towards  developing  more  effective  anisotropic  materials  with  larger  (and 
unequal)  strain  or  stress  constants  (d3i  or  cj?)  and  a  larger  modulus  of  elasticity. 

b)  The  finite  element  method  can  provide  an  accurate  solution  for  thin  and  thin-walled  structures 
microactuated  by  thin  piezoelectric  elements.  This  method  can  readily  account  for  the  presence  of  webs 
and  stiffeners,  either  active  or  passive,  to  study  new  configurations. 

Future  work  will  concentrate  on  using  this  finite  element  model  to  assess  aeroelastic  effects  on  thin  plates  in 
supersonic  flow  and  built-up  wings  in  transonic  flow.  Such  configurations  can  be  used  to  suppress  flutter,  control  shock 
wave  formation  and  to  reduce  drag.  In  addition,  this  finite  element  formulation  can  be  used  for  optimization  studies  to 
determine  effective  actuator  geometry  and  locations  and  to  help  to  understand  issues  related  to  integrating  these  actuators 
with  host  structures. 
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Abstract 

This  paper  examines  the  electromechanical  interaction 
between  a  thin  self-straining  piezoelectric  actuator  and 
a  simply  supported  host  plate  when  the  actuator  is 
placed  on  one  side  of  the  plate  and  its  objective  is  to 
create  a  large  bending  deflection.  The  purpose  of  the 
actuator  is  to  create  bending  deflections  to  control  local 
aerodynamic  pressures  and  resultant  forces  such  as  lift 
and  drag.  These  studies  show  that  there  is  a  trade-off 
between  the  additional  stiffening,  provided  by  actuator 
thickness  and  the  area  that  it  covers  on  the  plate,  and  the 
amount  of  force  and  moment  provided  by  the  actuator. 
A  Rayleigh-Ritz  analysis  shows  that  the  optimum  size, 
thickness  and  coverage,  of  the  actuator  with  respect  to 
the  host  panel  is  determined  by  panel  aspect  ratio,  and 
relative  elastic  moduli  The  strain  energy  content  of  the 
actuator/plate  combination  shows  that  the  best 
combinations  of  actuator  thickness  and  panel  coverage 
can  be  identified  by  plotting  strain  energy  against 
actuator  thickness  or  area.  With  an  aluminum  host  plate 
and  a  PZT  actuator,  the  best  rectangular  actuator  size  is 
about  0.6  the  thickness  of  the  plate  and  covers  about 
65%  of  the  host  plate. 

Imipductinn 

Aerodynamic  loads,  and  the  local  pressures  that 
create  these  loads,  depend  on  the  surface  shape  on 
which  they  act  Surface  panels  may  be  flat  or  curved 
and  are  designed  to  provide  aerodynamic  shape  and  to 
guarantee  structural  integrity.  Re-shaping  smooth 
aero/* true tural  surfaces  to  change  die  pressure 
distribution  is  done  by  bonding  or  otherwise  attaching 
thin  actuators  to  the  inside  surface*  of  structural  panels 
to  create  an  asymmetric  configuration  that  will  bend  on 
command. 

Thin  plate-like  or  lattice  reinforced  panels  with 
embedded  self-straining  actuators  such  ss  shape 
memory  materials  or  piezoelectric  materials  have  been 
proposed  for  aerodynamic  control  concepts  that  include 
actively  controlled  panels  to  reduce  transonic  drag1*2 
and  active  panel  elements  to  increase  supersonic  panel 
flutter  speed  3*4*5<&7.  For  transonic  drag  reduction,  the 
deformation  of  a  panel  on  the  upper  surface  of  a 
supercritical  airfoil  can  change  the  flow  field  and  shock 
wave  intensity  to  reduce  tkag  on  command. 
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Panel  flutter  suppression  with  piezoelectric 
actuators  and  shape  memory  alloy  actuators  is  unique  in 
that  no  vticulated  device  exists  to  do  the  same  ode  hi 
supersonic  flow,  dynamic  oscillations  can  be  reduced 
by  placing  thin  actuators  on  the  panel  surfaces  to 
change  the  frequencies  of  the  panels  on  command. 

One  serious  problem  with  active  panel  concepts  is 
the  difficulty  finding  a  design  combination  to  give  lane 
enough  panel  out-of-plane  deformations  to  create  the 
required  changes  in  aerodynamic  forces.  Without 
deflections  of  the  order  of  a  panel  thickness  (or  even 
more),  controlling  the  size  and  position  of  the 
aerodynamic  forces  is  marginal. 


A  desirable  actuator,  such  as  one  using  today's 
piezoelectric  materials,  can  not  create  significant 
bending  deformation  of  panels  unless  the  host 
psnd/actuator  combinabon  is  tailored  to  extract  every 
bit  of  electro-mechanical  efficiency  out  of  the 


configuration.  An  emphasis  on  efficiency  naturally 
leads  to  considerations  of  formal  optimization  that 


includes  a  design  objective  end  design  variables. 
However,  before  formal  optimization  can  proceed,  we 


must  select  our  design  variables  and  determine  the 
sensitivity  of  the  design  objective  to  these  design 


variables. 


This  paper  is  a  pre-optimization  study  that 
examines  the  interaction  between  actuator  self-straining 
ability,  bending  stiffness,  thickness  and  planform 
coverage  and  the  host  panel  bending  stiffness  and 
aspect  ratio.  The  purpose  of  the  actuator  is  to  produce 
large  bending  deflection.  The  intent  of  the  study  is  to 
identify  effective  panel/actuator  combinations  and 
understand  why  some  combinations  are  mare  effective 
than  others. 


Parameters  that  affect  panel  out-of-plane  deflection 
and  enter  into  the  optimization  process  include:  the 
thickness  of  the  actuator  compared  to  the  thickness  of 
the  panel;  the  surface  area  covered  by  the  actuator,  the 
position  of  the  actuator  on  the  panel  surface;  the 
boundary  conditions  at  the  edges  of  the  panel;  the 
aspect  ratio  of  the  panel  planform;  and,  the  material 
properties  of  the  actuator  compared  to  the  boat  panel. 

Analytical  work  described  in  this  paper  uses  a 
baseline  configuration  shown  in  Figure  1.  The  size  of 
this  panel  is  consistent  with  the  requirements  for  an 
Ktive  panel  that  might  be  placed  between  ribs  of  an 
active  wing.  The  boat  r'ate  material  is  aluminum,  with 
dimensions  of  12  in.  by  18  in.,  and  a  thickness  of  0.05 
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in.  A  thin  lead  zuxoaate  utanate  (PZT)  piezoelectric 
actuator  ts  attached  to  the  bottom  inner  surface  of  the 
host  plate.  Expansion  or  contraction  of  this  actuator 
will  cause  panel  extension  and  bending.  The  area  and 
thickness  of  the  actuator  are  parameters  for  the  study. 

Kim  A  Jones®  have  studied  a  similar  problem  to 
find  the  optimal  thickness  of  piezoactuators  to 
maximize  the  bending  moment  induced  by  rectangular 
actuators  surface-bonded  to  the  upper  and  lower 
surfaces  of  a  thin  flat  plate.  They  found  that  the  best  - 
actuator  thickness  is  approximately  half  the  thickness  of 
a  steel  host  substructure  and  a  quarter  of  (he  thickness 
for  an  aluminum  host  substructure. 


planform  view 


Figure  I  -  Pa  net/ actuator  geometry 

Rogers.  Liang.  A  Ji a9  studied  mid-plane  symmetric 
shape  memory  alloy  reinforced  plates,  using  a 
Rayleigh-Ritz  method  for  their  numerical  results.  They 
obtained  an  approximate  solution  to  (he  plate  bending 
problem,  free  vibration,  buckling,  and  acoustic 
transmission  loss. 

Crawley  A  Lazarus10  developed  a  consistent  plate 
model  with  embedded  actuator  stiffness  and  strain 
included  as  laminated  plate  layers.  Their  computations 
were  based  on  a  Rayleigh-Ritz  model.  Wang  A 
Rogers11  also  applied  laminate  plate  theory  to  a 
laminated  plate  containing  induced  strain  actuator 
"patches”  bonded  symmetrically  to  the  surface  or 
embedded  within  the  laminate.  The  thickress  and  size 
of  these  actuator  patches  are  relatively  small  compared 
to  those  of  each  lamina. 

This  study  differs  Grom  previous  studies  because 
the  actuator  is  placed  on  only  one  side  of  the  host  plate. 
However,  we  will  use  laminated  plate  theory  and  a 
Rayleigh-Ritz  solution  technique,  although  we  will 
check  its  accuracy  using  a  finite  element  analysis. 

U ns  ym  metric  host  plate/actuator  combinations  are 
different  than  symmetrical  combinations  in  at  least  two 


major  ways.  First  of  all,  nonsymmetry  of  the  material 
stiffness  through  the  thickness  of  the  plate/actuator 
combination  assures  that  there  will  be  coupling  between 
in -plane  induced  strain  and  panel  bending.  The  bending 
stiffness  of  the  plate/actuator  combination  will  always 
be  less  than  a  similar  combination  of  material  where  the 
total  actuator  thickness  is  the  same,  but  distributed 
symmetrically  about  the  mid-plane  *2 

There  is  a  second  major  difference  between  our 
configuration  and  symmetrical  panel  actuator 
combinations.  For  unsym metric  cross-ply  laminates, 
’large  deflection  effects’  can  occur  for  configuration 
loading  that  normally  would  be  regarded  as  producing 
deflections  in  the  small  plate  deflection  range.  Bending- 
extension  coupling  can  produce  a  stiffening  or  softening 
effect  depending  on  the  direction  of  the  deflection.1’ 
An  unsymmetric  cross-ply  laminate  in  cylindrical 
bending  has  different  apparent  bending  stiffnesses  in 
the  positive  and  negative  deflections.  As  a  result, 
positive  and  negative  loads  of  the  same  magnitude 

Eluce  different  magnitudes  of  deflection  and  linear 
maoon  theory  may  give  large  differences  when  the 
nonlinear  effects  ve  ignored.14 

Nonlinear  large  deflection  effects  will  most  likely 
increase  the  bending  deflection  found  from  linear 
theory.  Our  intent  is  to  provide  information  as  how  to 
maximize  panel  bending  deflections.  The  additional 
computational  effort  required  to  consider  these 
non  linearities  was  not  regarded  as  essential  to  these 
results  at  this  time,  so  nonlinear  effects  were  excluded 
from  this  study. 

Two  computational  methods  for  calculating 
bending  deflection  were  considered.  A  finite  element 
program  (NAS TRAN)  was  first  used  to  compute  the 
deformations  of  the  plate/actuator  combination.  In 
addition,  a  Rayleigh-Ritz  procedure  was  developed, 
beginning  with  a  strain  energy  expression  based  on 
laminated  plate  theory.  This  expression  and  the 
development  of  the  equations  of  static  equilibrium 
necessary  to  compute  bending  deflection  are  provided 
in  the  Appendix.  The  analysis  is  based  on  classical 
laminated  plate  theory.  The  energy  expression  is  useful 
as  a  guide  to  understanding  why  some  plate/actuator 
combinations  are  better  than  others. 

Finite  Element  Model 

The  panel  configuration  modeled  using  the 
NASTRAN  finite  element  program  consisted  of 
quadrilateral  plate  elements  with  membrane  and 
bending  stiffness.  Each  of  the  small  elements  has  the 
same  aspect  ratio  as  the  plate,  and  are  arranged  24  along 
each  side,  for  a  total  of  576  elements.  The  actuator 
elements  are  directly  attached  to  (he  host  plate  with  no 
intervening  bonding  layer. 

For  each  different  actuator /host  plate  configuration 
considered,  a  new  model  must  be  created  and  the  finite 
element  program  run.  While  the  model  computation 
time  is  small,  the  time  required  to  create  models  and  to 
interpret  data  was  considered  to  be  excessive,  given  the 
scope  of  our  study. 
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The  piezoelectric  strain  is  a  combination  of  the 
piezoelectric  constant  cf5l  multiplied  by  the  electric 
field  strength  £3.  NASTRAN  has  no  piezoelectric 
finite  element  capability,  so  an  equivalent  thermal 
element  was  used.  Thermal  and  piezoelectric  strains 
are  both  induced  strains  so  that  a  thermal  coefficient  of 
expansion  and  temperature  increase  can  be  assigned  to 
mimic  the  actuation  strain  of  the  piezoelectric  material. 


An  approximate  solution  for  the  plate  deflections  - 
can  be  obtained  using  a  Rayleigh-Ritz  method  based  on 
laminate  electro-mechanical  strain  energy  and  an 
assumed  displacement  field  for  inpiane  and  out-of-plane 
(bending)  deflections.  The  panel  strain  energy 
expression  in  the  Appendix  is  composed  of  four  basic 
types  of  terms  and  may  be  written  conceptually  as: 


The  reader  should  note  that  the  deflection  coefficients 
Am,  and  in  Eqns.  A  and  5  are  not  the  same  as  the 
laminate  stiffness  coefficients  in  Eqn.  1 . 

When  either  of  the  assumed  displacement  expressions 
are  substituted  into  the  strain  energy  expression,  the 
Principle  of  Virtual  Work  can  be  used  so  that  the  strain 
energy  is  minimized  with  respect  to  the  displacement 
coefficients.  A  set  of  linear  equations  of  static 
equilibrium  results. 

When  the  polynomial  deflection  is  used,  this 
energy  minimization  results  in  only  three  equilibrium 
conditions  for  the  constants  A.  B  and  C.  given  as 
follows 


<J) 


U  =  i nplane  stiffness  +  bending  stiffness 
+ uiplane  /  bending  coupling  +  induced  strain 

The  stiffnesses  involving  shear-extension  coupling 
vanish  if  the  laminate  is  a  symmetric  laminate  with 
isotropic  or  specially  orthotropic  layers,  that  is 

=  Ajj  =  Big  *  *  0)6  *  ^26  a  0  (1) 

The  assumed  displacements  must  satisfy  the 
simply-supported  plate  edge  boundary  conditions.  Two 
approximate  solutions  were  used.  The  first  is  a  single 
polynomial  term  (a  true  Rayleigh  solution),  while  the 
other  is  a  mare  general  series  solution  (a  Ritz  solution). 

When  the  actuator  is  restricted  to  be  at  the  plate 
center,  a  single  set  of  three  polynomial  assumed 
displacements,  with  their  origin  at  the  center  of  the 
plate,  is 


This  set  of  three  simultaneous  equation  sets  must  be 
solved  for  the  deflection  coefficient  sets  a.  B.  and  c 
These  coefficients  are  then  substituted  into  the 
displacement  equations  to  find  the  plate  deflection  at 
the  center  of  the  plate. 

If  only  bending  of  the  plate  is  considered,  and  the 
inplane  energy  and  inplane-bending  coupling  energy  is 
ignored,  the  strain  energy  expression  contains  only 
tending  stiffness  matrix  terms  and  the  induced 
piezoelectric  moments.  To  solve  for  the  deflections, 
only  the  third  of  the  polynomial  equations  is  subarituted 
into  the  energy  expression.  While  this  reduces  the 
workload,  it  leads  to  errors.  The  coupling  terms  are 
very  important. 

There  are  differences  between  retaining  the  inplane 
stiffness  and  inplane/bending  stiffness  coupling  terms 
and  ignoring  item.  Figure  2  compares  the  results  of  the 
uncoupled  tending  analysis  and  the  bending-extension 
coupling  analysis  using  the  simple  polynomial 
expression  to  the  finite  element  solution.  All  of  the 
were  run  far  an  artuatnr  10  plate  thickness  ratio  of 

0.6. 


where  u.  v  and  w  are  the  inplane  displacements  in  the 
x.y  and  tending  directions,  respectively. 

A  more  general  assumed  displacement  set  is  ^ 

“o  ”  X  X  (4) 

«*>1  »"l 

%  =  X  X  (5) 

«■!  *■! 


Comparison  ot  Rayteigh-Altz  to  NASTRAN 


Actator  Area /Pit*  Atm 


Figure  2  •  Plate  ceatcr  bending  deflection  coni  pa  ted 
by  Rita  Method  with  polynomial  approximation. 
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It  can  be  seen  that  the  uncoupled  bending  solution 
over  predicts  the  deflection  for  the  entire  range  of 
actuator  area  ratios,  showing  that  it  is  important  to 
include  the  inplane/bending  coupling  terms.  The 
bending-extension  coupling  solution  is  remarkably 
close  to  the  NASTRAN  solution  for  these  simple 
displacement  functions.  However,  the  actuator  ij 
restricted  to  being  centered  on  the  plate. 

Because  we  want  to  assess  the  effects  of  changing 
the  actuator  location  on  the  plate,  and  to  get  a  more 
accurate  solution,  the  trigonometric  series  displacement 
functions  were  used  for  tall  other  studies.  Substituting  - 
them  into  the  strain  energy  expression  in  the  Appendix 
and  minimizing  with  respect  to  the  undetermined 
displacement  coefficients  results  in  the  general 
conditions 


dU 


0 


(8) 


The  result  is  a  set  of  3  M  x  N  simultaneous  equations 
which  must  be  solved  for  the  series  coefficients  A^,, 
and  C^. 


The  form  of  the  equations  for  the  displacement 
coefficients  is 


(9) 

where  is  a  (3  M  x N)X(3 Af x N)  matrix  of 

stiffness  terms,  found  from  integrating  the  assumed 
displacement  series  strains  into  the  energy  expression. 

{Am}  ua  3  MxN  vector  consisting  of  the  vectors  of 

Am,  0M,  and  C^,  coefficients,  (fiy}  is  a  3A#xAf 

load  vector  whose  terms  represent  the  induced 
piezoelectric  forces  and  moments.  These  latter  terms 
can  be  thought  of  as  the  modal  farces  and  moments  for 
each  of  the  assumed  displacement  terms. 

The  stiffness  matrix  [ATy,**]  is  inverted  to  solve  for 

the  displacement  coefficients..  This  matrix  becomes 
large  as  more  terms  are  added  to  the  series  solution. 
For  our  results.  M  =  N  - 11  was  used.  This  results  in 
a  stiffness  coefficient  matrix  that  is  of  order  363x363. 
Figure  3  shows  the  agreement  between  this  series 
solution  and  the  finite  element  results. 

FffrvM  n i  mm  ifajdaca 

Actuator  thickness  and  area  far  the  baseline  plate  in 
Figure  1  were  varied  to  find  an  actuator/riuminum  host 
plate  combination  that  produced  the  moat  panel  center 
bending  deflection. 

For  these  studies,  the  actuator  was  located  at  the 
center  of  the  plate  planform  and  attached  to  its 
underside.  The  maximum  bending  deflection  occurred 
at  the  plate  center.  A  range  of  actuator  thickness  to 
plate  thickness  ratios  between  0.1  to  1.0  was  examined. 
For  each  thickness  ratio,  a  range  of  actuator  planform 
coverage  from  3%  to  100%  of  the  plate  area  was 


examined.  From  these  results,  the  combination  of 
actuator  thickness  ratio  and  area  redo  that  produced  the 
largest  deflection  could  be  determined. 
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Ftgnra  3  •  Ptaia  cur  bending  defection  arias 
trigonometric  Ma  apprnrimatton  with  thkknaaa 
ratio  of  04  •  comparison  with  NASTRAN  (M>N*I1) 

Figure  4  shows  the  results  of  a  Rayleigh-Ritz 
analysis  to  compute  the  panel  center  deflection 
produced  by  an  artuator  with  an  actuator/host  plate 
thickness  ratio  between  0.1  and  0.6  for  actuator/host 
plate  area  ratio  between  0. 1  and  1.  This  figure  indicates 
that  an  actuator/plaie  thickness  ratio  of  0.6 produces  the 
largest  j»el  deflection  when  the  actuator/piate  area 


Comparison  of  Sartos  Assult  to  NASTRAN 


The  curvet  in  Figure  4  show  that  a  slightly  smaller 
area  ratio  than  0.6  will  produce  nearly  as  much  panel 
center  deflection.  Using  a  smaller  actuator  would 
produce  a  weight  savinp  with  little  degradation  in 
performance.  An  actualor/plate  thickness  ratio  of  0.S 
also  produces  its  greatest  deflection  at  65%  area 
coverage.  Note  that  the  deflection  produced  with  an 
actuator  0J  as  thick  as  the  plate  is  not  much  different 
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than  that  produced  by  an  actuator  whose  thickness  ratio  changing  the  actuator  aspect  ratio  on  the  panel  center 
is  0.6.  bending  deflection  was  examined  using  our  model. 


When  the  actuator  has  a  thickness  ratio  of  0.1.  a 
maximum  in  the  curve  of  deflection  vs.  actuator 
coverage  occurs  at  about  804  coverage.  However,  the 
deflection  does  not  change  much  between  50-100% 
coverage.  As  the  actuator  thickness  increases,  the 
maximum  value  in  the  panel  center  deflection  curve 
becomes  more  pronounced  and  shifts  toward  a  smaller 
actuator  area  coverage 

Figure  S  shows  the  results  of  analyses  for 
aciuator/plate  thickness  ratios  between  0.6  and  1  for  the 
enure  area  ratios  between  0.1  and  1.  The  thickness  ratio 
of  0.6  still  produces  the  most  deflection  at  area  ratios 
greater  than  50%.  However,  thickness  ratios  of  0.7  and 
0.8  produce  larger  deflections  at  low  actuator  area 
coverage. 


Figure  5-  Plate  center  beoding  deflection  v*. 
actuator  area  coverage  for  thkkaeas  ratios  betweea 
0.6  and  lit  (iacrements  of  0.1);  aspect  ratio  1-5. 


Effects  of  panel  aspect  ratio 
on  oodmal  actuator  size 


For  the  square  panel,  the  best  actuator  has  an 
actuator  thickness  ratio  of  0.6  and  an  area  coverage  of 
65%.  Any  actuator  aspect  ratio  less  than  or  greater  than 
1  creates  a  smaller  center  deflection  than  the  initial 
configuration,  although  the  differences  are  small. 

When  the  panel  aspect  ratio  is  l  .5,  the  best  actuator 
also  has  a  thickness  ratio  of  0.6  and  coven  an  area  of 
65%  of  the  panel.  When  the  actuator  aspect  ratio  is 
increased  so  the  actuator  spans  the  long  dimension  of 
the  panel,  the  greatest  decrease  in  deflection  is  about 
5.7%. 
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Figure  6  -  The  effect  of  actuator  coverage  oa  center 
deflection  of  three  different  pane  la;  constant 
actuator  and  panel  area  with  actna tor/plate 
thickness  ratio  0.6;  aspect  ration  1, 1-5  and  3. 

For  a  panel  with  an  aspect  ratio  of  3.  the  best 
actuator  has  a  thickness  ratio  of  0.6  and  an  area 
coverage  of  70%.  For  this  case,  an  actuator  with  an 
aspect  ratio  smaller  than  the  panel  produces  the  greatest 
deflection,  although  the  increase  is  only  0.5%. 


Figure  6  shows  the  maximum  panel  center  bending 
deflection  obtained  at  an  actuator/host  plate  thickness 
ratio  of  0.6  for  several  different  panel  aspect  ratios. 
Although  the  aspect  ratio  can  be  changed,  the  total  plate 
area  and  actuator  area  remain  fixed  to  be  the  same  as  a 
panel  with  dimensions  18  in.x  18  in.  The  square  panel 
with  an  aspect  ratio  of  1  has  the  largest  center 
deflection.  The  panel  with  an  aspect  ratio  of  1.5  has  a 
deflection  smaller  than  that  of  the  square  panel,  but  the 
maximum  bending  deflection  still  occurs  when  the 
actuator  coverage  is  about  65%  of  the  plate. 

The  panel  with  a  planform  aspect  ratio  of  3  has  a 
maximum  center  deflection  less  thsn  half  of  the  square 
panel,  even  though  the  panel  areas  are  the  same.  In  this 
case  the  simple  supports  are  so  close  to  the  panel  center 
when  the  aspect  ratio  is  3  that  the  stiffness  of  the  panel 
is  increased  with  respect  to  center  deflection. 

For  all  of  the  previous  studies,  the  actuator  area  and 
aspect  ratio  were  restrained  to  be  the  same  as  the  plate 
to  which  it  was  mounted,  although  die  actuator  aspect 
ratio  does  not  need  to  be  the  same.  The  effect  of 


Nnn-fflunilY  Inraml  acnamn 

For  all  previous  results,  the  actuator  was  centered 
on  the  panel  planform  and  the  maximum  deflection 
occurred  in  the  center  of  the  panel.  For  a  non-centered 
actuator,  the  maximum  deflection  will  not  occur  in  the 
center  of  thepsnel,  so  a  deflection  distribution  must  be 
examined.  This  centerline  is  located  in  the  x -direction, 
parallel  to  the  long  edge  of  the  plate. 

The  baseline  panel  with  an  aspect  ratio  of  1.5  was 
used  for  a  study  of  the  effects  of  non  centrally  located 
actuators.  The  actuator  has  the  same  aspect  ratio  as  the 
panel  itself.  Two  different  actuator  sizes  were  used; 
they  covered  6.25%  and  25%  of  the  plate  ares.  The 
actuator  was  subjected  to  its  maximum  electric  field 
throughout  and  located  at  one  of  four  different 
locations  along  the  centerline.  In  each  case,  the 
deflection  profile  was  calculated  and  compared. 

The  smallest  (6.25%  area  coverage)  actuator  was 
placed  in  the  center  of  thepeneL,  then  shifted  toward  the 
left  edge  of  the  panel.  Figure  7  shows  the  deflection 
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profile  at  the  centerline  of  the  panel  in  each  case  (.see 
inset  for  positions).  The  actuator  edge  can  be  moved 
close  to  the  edge  of  the  panel  without  seeing  a  large 
decrease  in  the  maximum  deflection.  However,  once  an 
actuator  edge  is  at  or  near  the  panel  edge,  the  peak 
deflection  decreases  noticeably.  Since  the  shape  and 
amplitude  of  the  bending  deflection  on  the  panel  can  be 
changed  by  relocating  the  actuator  this  provides  a  way 
to  tailor  the  aerodynamic  shape. 


Panel  with  Actuator  on  6.25%  of  Area 
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Figure  7  •  Bending  deflectioa  distribution  for  a 
panel  with  an  actuator  covering  of  panel. 

Figure  8  shows  the  centerline  deflection  profile  far 
the  actuator  covering  25%  of  the  panel  area  and  how 
this  deflection  changes  as  the  actuator  is  repositioned 
along  the  panel  centerline.  The  maximum  deflection 
here  is  about  twice  that  of  the  smaller  6.25%  coverage 
actuator.  The  peak  deflection  when  the  actuator  edge  is 
at  the  edge  of  the  panel  is  about  16%  less  than  when  the 
actuator  is  centered. 
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Figure  8  -  Bending  deflection  distributions  for  a 
panel  with  an  actuator  covering  25%  of  panel  area. 

Strain  energy  denary  for  optimal  actuator  size 

We  have  observed  that  actuator  thickness  and  panel 
area  coverage  affect  the  size  of  the  panel  bending 


Panel  with  Actuator  on  25%  of  Area 
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deflection.  When  we  try  to  create  bending  deflection 
with  this  type  of  actuator,  two  effects  are  in  conflict 
with  each  other.  Fust  of  all.  as  actuator  thickness  (or 
area)  increases  the  actuator  is  able  to  create  large  strains 
to  induce  inplane  forces  and  bending  moments.  On  the 
other  hand,  as  the  actuator  becomes  thicker  and  coven 
more  area,  the  panel  stiffness  also  increases  so  that  the 
plate  is  more  difficult  to  bend. 

When  the  actuator  is  small  compared  to  the  plate, 
its  attempt  to  expand  is  easily  thwarted  by  a  relatively 
massive  plate.  When  a  voltage  is  applied  to  the 
actuator,  most  of  the  strain  energy  comes  from  the 
actuator  being  held  relatively  fixed  and  expanding  only 
inplane  where  the  inpline  stiffness  is  larger  than  the 
sc tua lor  bending  stiffness.  Very  little  strain  energy 
goes  into  bending. 

On  the  other  hand,  as  the  actuator  grows  in  size,  the 
rate  of  change  of  induced  bending  moment  with  respect 
to  actuator  thickness  is  relatively  large  compared  to  the 
rate  of  change  of  panel  bending  stiffness.  As  a  result,  it 
is  better  to  increase  the  size  of  the  actuator  and  more 
energy  goes  into  bending.  Bending  stiffness  is  much 
less  than  inplane  stiffness  so  the  rale  of  change  of  strain 
energy  with  respect  to  changes  in  actuator  thinking 
becomes  less. 

As  the  actuator  thickness  increases,  the  bending 
stiffness  increases  more  rapidly  than  the  induced 
moment.  At  some  point,  the  rate  of  change  of  induced 
bending  moment  will  equal  the  rate  of  change  in 
bending  stiffness.  This  is  the  best  actuator  for  bending 
deflection.  Any  increase  in  actuator  size  will  increase 
bending  stiffness  more  than  induced  bending  moment 
and  will  be  counter-productive.  The  actuator  becomes 
non-optimaL 

The  dependence  of  panel  strain  energy  on  actuator 
thickness  can  be  plotted  for  a  panel  with  two  parallel 
edges  free  while  the  other  two  edges  are  simply 
supported.  This  is  the  case  called  cylindrical  bending. 
When  the  actuator  coven  the  entire  plate  area,  the 
resulting  deflection  when  the  actuator  is  operated  at  its 
full  power  is  a  parabola  with  constant  curvature.  The 
expression  for  strain  energy  of  this  special  plate  is  as 
follows: 


Consider  the  case  of  a  constant  thickness  aluminum 
plate  with  a  PZT  actuator  covering  one  100%  of  the 
plate.  The  strain  energy  (U  divided  by  the  area  of  the 
plate)  is  plotted  in  Figure  9  against  actuator  thickness 
ratio.  Also  plotted  is  the  curvature  of  the  plate  vs. 
thickness  ratio. 

The  maximum  bending  deflection  occurs  when  the 
plate  curvature  is  a  maximum  and  when  the  actuator 
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thickness  ratio  is  0.5.  The  strain  energy  density  has  an 
inflection  point  at  this  thickness  ratio. 


Figure  9  -  Btsa  strata  energy  density  aad 
curvature  vs.  actuator/plate  thickness  ratio.  Full 
coverage  actuator;  cylindrical  bending. 


The  presence  of  an  inflection  point  is  important 
because  it  suggests  an  optimality  criterion  to  select  the 
actuator.  The  ’too  small”  actuator  stores  a  large  part  of 
the  mechanical  energy  because  the  panel  barely 
deforms.  When  the  actuator  reaches  a  'critical  size*  the 
panel  bends  a  great  deal,  but  stores  energy  in  a  more 
'flexible'  mode.  As  actuator  sue  increases  further,  the 
dominant  energy  storage  mode  again  becomes 
extensional  and  the  slope  of  the  strain  energy  density 
curve  turns  upward  and  increases  rapidly. 

The  favorable  bending  energy  storage  mode  is 
identified  by  an  inflection  point  or  ’flat  spot’  in  the 
strain  energy  density  vs.  actuator  thickness  curve. 


Plata  Strain  Energy  for  3  Area  Ratiea 


aspect  ratio  is  1.5. 


The  strain  energy  for  the  plate  can  also  be  plotted 
versus  thickness  ratio.  Figure  10  shows  the  strain 
energy  plotted  vs.  thickness  ratio  for  actuators  covering 
609b,  65%,  and  70%  of  the  plate  area  for  an  aluminum 
plate  with  an  aspect  ratio  of  I.S.  In  the  cylindrical 
bending  case,  it  was  possible  to  identify  the  best 
thickness  ratio  by  looking  for  an  inflection  point  in  the 
strain  energy  plot  However,  this  inflection  point  is  not 
so  evident  far  the  simply  supported  panel. 

The  curves  in  Figure  10  become  nearly  linear  for 
larger  thickness  ratios.  Figure  10  indicates  that,  as  long 
as  the  slope  of  the  strain  energy  curve  is  changing  with 
increases  in  thickness  ratios,  adding  thickness  to  the 
actuator  will  be  beneficial.  But  once  the  curve  becomes 
linear,  adding  actuator  thickness  will  not  produce  more 
bending  deflection  in  the  panel. 

Conchrama 

A  panel  with  a  self-straining  actuator  mounted  on 
one  side  was  studied  to  find  the  features  of  the  actuator 
that  produces  the  largest  deflection  of  simply  supported 
rectangular  panels.  A  Rayleigh-Ritz  model  was 
developed  to  compute  inplane  and  betiding  deflections 
of  a  plate  with  an  actuator  covering  only  pan  of  the 
area. 


It  was  shown  that  it  is  important  to  include  the 
coupling  terms  and  the  in-plane  actuator  forces  to 
accurately  model  the  problem.  A  simple  polynomial 
assumed  displacement  field  provided  good  agreement 
with  the  finite  element  results  for  centrally  located 
actuators  with  a  large  area  coverage. 

To  match  the  finite  element  results  far  smaller 
actuator  areas  and  to  provide  more  accurate  results 
overall,  a  Rayleigh-Ritz  trigonometric  series  expansion 
was  developed.  This  model  allows  the  actuator  to  be 
placed  at  any  location  on  the  panel.  It  was  found  that  a 
series  with  M  =  S  =  \\  terms  far  each  of  the  u,  v.and 
w  displacements [yielded  excellent  results,  very  close  to 
those  of  the  NASTRAN  program. 


For  aluminum  panels  with  aspect  ratios  between 
1.0  and  15,  the  best  actuator  has  a  thickness  ratio  of  0.6 
and  coven  65%  of  the  panel  area.  For  a  panel  with  an 
aspect  ratio  of  3.0,  the  best  actuator  has  a  thickness  ratio 
of  0.6  and  coven  70%  of  the  area.  For  the  panels  with 
aspect  ratios  of  1.0  and  15.  the  actuator  with  the  same 
aspect  ratio  as  the  panel  produced  the  largest  deflection. 
For  the  panel  of  aspect  ratio  3.0.  an  actuator  with  an 
aspect  ratio  slightly  smaller  than  that  of  the  panel 
produced  the  most  deflection. 

The  effect  of  actuator  location  on  the  deflection 
profile  of  a  simply  supporied  panel  was  examined.  An 
actuator  with  25%  area  coverage  or  less  produced  the 
most  change  in  the  bending  deflection  shape. 

For  a  plate  with  cylindrical  bending,  when  the 
strain  energy  is  plotted  against  thickness  ratio,  the  curve 
shows  an  reflection  point  at  the  best  actuator  thickness 
for  that  configuration.  Although  the  best  actuator 
thickness  was  obviously  at  the  reflection  point  in  the 
strain  energy  curve  for  the  cylindrical  bending  case,  the 
best  actuator  thickness  for  the  general  plate  was  less 
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best  actuator  thickness  for  the  general  plate  was  less 
obvious. 

The  plot  of  stnun  energy  vs.  actuator  thickness 
ratio  becomes  nearly  linear  at  the  point  corresponding 
to  the  best  actuator  thickness  ratio.  This  indicated  that 
increasing  the  thickness  of  the  actuator  was  beneficial 
until  the  slope  of  the  strain  energy  curve  reached  its 
smallest  value.  Adding  more  thickness  to  the  actuator 
beyond  this  point  increases  stiffness  more  than  it 
increased  the  applied  moments,  and  does  not  increase 
the  deflection  dial  is  obtained  in  the  panel. 

For  this  study,  actuator  thickness,  area,  and  aspect 
ratio  were  varied  to  find  which  combination  produced 
the  largest  deflection  in  a  given  panel.  More  precise 
results  for  the  actuator  characteristics  could  be 
calculated  if  an  optimization  scheme  were  used  in 
conjunction  with  the  model  already  developed.  Since 
only  a  finite  number  of  combinations  were  tried  in  this 
study,  the  best  actuator  was  found  within  the  limits  of 
the  study.  Also,  the  optimal  actuator  could  be  found  by 
taking  into  consideration  the  weight  that  is  added  for  a 
larger  actuator.  A  larger,  thicker  actuator  does  not 
always  produce  significantly  more  deflection  than  a 
smaller,  lighter  actuator. 

This  study  was  also  limited  to  rectangular  actuators 
placed  at  any  location  on  the  panel.  Actuators  of  stupes 
other  than  rectangular  should  also  be  considered. 
Rectangular  actuators  produce  high  stresses  at  the  edges 
and  corners17.  The  comer  stresses  might  be  avoided  if 
the  actuator  shape  were  changed.  The  deflection 
produced  by  elliptical  actuators  or  other  actuator  shapes 
should  be  examined. 

The  plate  deflections  obtained  here  used  the  PZT 
actuators  to  their  fullest  extent,  applying  voltages  to  the 
actuators  that  are  right  at  their  capabilities.  Much  needs 
to  be  done  to  improve  the  efficiency  of  these  materials. 
On  the  other  hand,  we  do  believe  that  there  is  promise 
for  future  applications  to  high  speed  flow  control. 
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Appendix  ■  Strain  Energy  Exnre«rinn 


a  plate  length 

\  extensions!  stiffness  matrix 

A,  A^,  assumed  series  solution  constants 

b  plate  width 

flv  coupling  stiffness  matrix 

fl.fl^i  assumed  series  solution  constants 
C.C^  assumed  series  solution  constants 
d]t  piezoelectric  constant 
Ov  bending  stiffness  matrix 
£j  electric  field  strength 


N*.'N*<Nw*  piezoelectric  inplane  forces 
piezoelectric  moments 

Gy*’  reduced  lamina  stiffness  far  kih  layer 

Tt  actuator  thickness 

T,  host  plate  thickness 

u.  inplane  defkctian  in  x  direction 

U  panel  min  energy 

v,  inplane  deflection  m  y  direction 

w  bending  deflection 

A(*>  piezoelectric  strain  of  kih  lamina 


The  strain  energy  expression,  U ,  is  used  to  oompuse  the  mechanical  strain  far  a  plate/actuaior  combination  in 
plane  stress  when  the  piezoelectric  actuator  creates  induced  strain.  The  sram  is  the  difference  between 

the  total  strain  and  the  expansion  strain,  and  is  the  only  strain  that  cream  sores 16.  There  is  no  sores  if  the 

material  is  allowed  to  expand  finely.  The  induced  piezoelectric  strains  are  analogous  to  thermal  strains. 

To  compute  strain  energy,  a  reference  surface  at  the  mid-plane  of  the  host  plate  was  chosen.  This  allows  the 
reference  surface  to  remain  fixed  when  the  actuator  thickness  changes  thus  reducing  computational  complexity.  The 
last  three  terms  of  the  equation  are  independent  of  the  displacements  it,  v,  and  w.  The  integral  involving  these 
terms  will  vanish  under  the  first  variation  of  the  strain  energy  and  the  terms  will  not  enter  the  equations  when  the 
Rayleigh-Ritz  method  is  vptied.  For  this  analysis  a  rectangular  actuator  is  attached  to  a  plate  at  re  sbitrary  position. 
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